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Positions  and  proper  motions  for  113  stars  in  the  Orion 
Nebula  are  derived  from  measurements  on  92  plates  taken 
with  a  variety  of  telescopes,  with  an  average  epoch  difference 
of  50  years.   The  positions  and  motions  are  measured  with  re- 
spect to  a  set  of  some  40  reference  stars;  they  are  therefore 
instrumentally  absolute  and  unaffected  by  changes  in  plate 
scale  over  time.   In  order  to  incorporate  a  sufficient  num- 
ber of  reference  positions  into  the  solution  and  to  obtain 
the  greatest  possible  accuracy  from  the  data,  the  reductions 
were  carried  out  by  the  plate  overlap  method.   A  new  version 
of  the  method,  particularly  economical  for  the  astrometry  of 
star  clusters,  was  employed. 

It  is  found  that  1)  the  velocity  dispersion  of  the  stars 
in  the  Trapezium  Cluster  (within  6'  of  the  Trapezium)  is  less 
than  2%  km/sec,  and  2)  the  Orion  Nebula  Cluster  (radius  20') 
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is  contracting  at  a  linear  rate  of  (V009/yr) /degree ,  cor- 
responding to  10  km/sec  at  the  4  parsec  radius  of  the 
radio  molecular  cloud.   These  results  imply  that  the  Orion 
Nebula  is  free  from  large  turbulent  motion,  is  gravitationally 
contracting,  and  that  star  formation  there  has  probably  been 
proceeding  continuously  since  it  began,  with  the  greater 
period  of  star  formation  yet  to  come.   These  conclusions 
in  turn  are  in  agreement  with  the  large  (-10^  solar  masses) 
mass  now  assigned  to  the  Orion  Nebula. 
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SECTION  I 
INTRODUCTION 


Background 
The  Orion  Nebula  has  been  an  object  of  special  fascina- 
tion to  astronomers  from  its  discovery  three  centuries  ago 
up  to  the  present.   It  is  now  recognized  as  a  complex  in- 
volving a  ("the")  trapezium,  a  star  cluster,  emission 
nebulosity,  infra-red  sources  and  radio  molecular  clouds. 
The  coincidence  of  these  components  in  a  single,  relatively 
nearby  (~500  pes)  object  potentially  makes  it  a  key  to 
understanding  a  number  of  basic  as  trophy sical  processes. 
At  the  same  time,  there  are  some  perplexing  problems  which 
stand  in  the  way  of  arriving  at  a  unified  model  of  the 
object . 

Let  us  examine  some  of  the  special  features  of  the 
Orion  Nebula  and  the  questions  relating  to  them. 

1.  The  Trapezium  itself  is  an  interesting  object,  the 
prototype  of  a  small  class  of  systems  ("Trapezium  system") 
first  recognized  by  Ambartsumian  (1955).  These  systems  are 
intermediate  between  binaries  and  true  clusters.  Even  when 
a  trapezium  is  considered  as  part  of  a  surrouiiding  cluster, 
it  is  clear  we  are  dealing  with  a  distinct  morphological 
group  since  the  trapezium  always  stands  out  as  a  small  (4-8 
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members) ,  central  grouping  at  least  five  magnitudes  brighter 
than  the  next  brighter  cluster  members,  in  contrast  to  ordi- 
nary open  clusters.   Trapezia  are  known  to  be  young  because 
they  occur  only  among  0  and  B  stars.   They  would  then  be 
disrupted  within  about  2  to  3.10   years,  yet  astrometric 
studies  fail  to  decide  whether  the  Orion  Trapezium  has 
positive  or  negative  energy. 

According  to  current  thinking  on  the  formation  of 
trapezia,  they  should  consist  exclusively  of  0  and  B  stars, 
yet  the  Orion  Trapezium  has  two  "comites"  (O'E  and  e'F) 
both  three  magnitudes  fainter  than  the  faintest  of  the  four 

principal  members. 

2.   The  Trapezium  is  surrounded  by  a  faint  cluster,  the 
"Trapezium  Cluster"  (or  "Orion  Nebula  Cluster"),  discovered 
only  as  late  as  1931  (Trumpler,  1931;  Baade  &  Minkowski, 
1931).   Its  brightest  members  are  ll~th  mag.  and  it  is 
usually  observed  by  the  surrounding  nebulosity.   The  clus- 
ter shows  well  in  the  plate  reproduced  in  Figure  1;  on  the 
original  plate  the  nebulosity  was  deliberately  suppressed. 

There  are  widely  differing  values  given  for  the  dia- 
meter of  this  cluster.   Baade  and  Minkowski  judged  its  dia- 
meter as  3'  to  5',  encompassing  some  100  stars.   Sharpless 
(1966),  referring  to  the  Baade  and  Minkowski  discovery, 
misquotes  the  diameter  as  "10'";  actually,  the  10'  figure 
appears  in  another  context  in  another  paper  by  the  same 
authors  immediately  preceding  (in  the  journal)  the  cluster 


FIGURE  1 
THE  ORION  NEBULA  SHOWING  THE  TRAPEZIUM  CLUSTER 


Two-fold  enlargement  of  plate  taken,  on  103a-G  emulsion 
with  3-69  filter,  exposed  for  12  min.  at  University  of 
South  Florida  Observatory.   A  circle  1*  in  radius  has  been 
drawn  around  the  star  e'C.   The  faintest  stars  visible  are 
14-th  magnitude. 


discovery  paper.   Trumpler,  on  the  other  hand,  though  count- 
ing the  stars  only  in  A'  diameter  circle,  actually  gives 
its  diameter  as  15'  -  partly  on  morphological  grounds. 
Gradually,  this  larger  value  has  crept  into  the  literature. 
Strand  (1958)  takes  the  diameter  as  30'. 

Turning  to  the  evidence,  one  can  see  in  Fig.  1  a 
definite  concentration  of  stars  within  a  circle  of  2' 
diameter  immediately  surrounding  the  Trapezium.   We  shall 
term  this  group  the  "Trapezium  cluster."   It  is  clearly 
much  smaller  than  a  typical  open  cluster,  and  morphologi- 
cally quite  different;  perhaps  it  might  be  better  named  the 
"Extended  Trapezium."   A  star  count  reveals  a  further  con- 
centration of  stars  in  an  area  of  about  5'  radius,  extending 
north  and  also  south-east  toward  the  group  6^.   It  is  to 
be  noted  that  heavy  absorption  obscures  the  region  south- 
east of  the  Trapezium.   In  the  absence  of  the  dust  we  should 
expect  to  see  cluster  members  in  that  quadrant  as  well.   (It 
is  to  be  expected  too  that  many  of  the  few  stars  visible  in 
that  area,  extending  for  at  least  20',  are  foreground 
objects.)   Accordingly,  we  assign  a  diameter  of  12'  to  this 
grouping,  calling  it  the  "inner  Orion  Nebula  Cluster." 
Finally,  in  accord  with  current  usage,  we  apply  the  name 
"Orion  Nebula  Cluster"  to  the  whole  region  of  >30'  in  diameter 
centered  on  the  Trapezium,  realizing  that  the  stars  in  this 
larger  region  may  not  comprise  as  well  demarcated  a  group. 


3.  Separate  theoretical  arguments  by  Menon  (1963,  Kahn 
and  Menon  1961)  and  by  Vandervoort  (1963)  ,  based  on  the  mass 
motions  and  ionization  of  the  inner  region  of  the  nebula, 
both  lead  to  an  age  of  only  Z-IO'*  years  (the  "ultra-short" 
age)  for  the  Trapezium  stars.   This  age  refers  to  the  time 
since  the  stars  have  been  capable  of  ionizing  the  hydrogen. 

4.  The  proper  motions  of  the  cluster  stars  can  in 
principle  yield  the  cluster's  age,  if  the  motions  show  a 
radial  expansion.   Then  the  kinematic  age  is  the  reciprocal 
of  the  expansion  constant.   (Strictly,  this  age  refers  only 
to  the  time  since  greatest  concentration  of  the  cluster  and 
does  not  preclude  earlier  contraction  and  star  formation.) 

Some  astrometric  studies  (Parenago  1953;  Franz  unpub- 
lished, quoted  by  Sharpless  1966)  have  indeed  confirmed  the 
ultra-short  age  of  the  Trapezium  and  this  figure  (lO'*  yrs) 
has  become  widely  quoted  in  the  literature.   But  other  astro- 
metric  studies  (Strand  1958,  Duboshin  et  al  1971)  result  in 
the  much  longer  (though  still  very  young)  kinematic  age  of 
2«10^  years;  while  yet  another  (Meurers  1963)  shows  no  evi- 
dence of  contraction,  and  still  others  (Cannell  unpublished, 
Vaerewyck  1972)  report  a  slight  contraction.   Adding  to  the 
confusion,  some  (Meurers,  Strand)  of  these  studies  show  a 
rotation  of  the  cluster.   The  compendious  work  of  Parenago 
(1954)  has  not  yet  been  analyzed  for  evidence  of  cluster  ex- 
pansion.  Parenago's  material  is  extremely  heterogeneous: 
significantly,  the  stars  for  which  he  quotes  the  smallest 
proper  motion  errors  also  have  the  smallest  dispersions  in 
proper  motions. 


5.   Color-magnitude  diagrams  for  the  Orion  Nebula 

1       2 
Cluster  (including  6   and  6  )  and  surrounding  regions  as 

given  by  several  investigators,  while  disagreeing  among 
themselves,  agree  in  pointing  to  a  young  age  for  the  stars 
in  Orion.   Some  of  the  discrepancy  can  be  attributed  to  the 
different  sizes  of  the  area  surveyed,  while  some  is  also 
due  to  the  differing  reddening  corrections  applied.   Walker 
(1969)  presents  a  particularly  careful  study.   His  photom- 
etry  shows  the  zero-age  main  sequence  down  to  an  intrinsic 
(B-V)  =  -.07,  corresponding  to  a  contraction  time  of 
~2  2/3  •  10   years.   At  fainter  magnitudes,  the  scatter  of 
points  above  the  main  sequence  would  seem  to  suggest  non- 
coevality  of  star  formation,  but  some  of  the  scatter  is  no 
doubt  due  to  the  patchy  distribution  of  the  absorption. 

Significantly,  no  photometry  has  been  carried  out  on 
the  stars  of  the  Trapezium  Cluster  (except,  of  course,  for 
the  Trapezium  itself) ,  due  to  the  difficulty  of  observing 
through  the  nebulosity.   Its  color-magnitude  diagram  might 
well  differ  from  that  of  the  larger  "Nebula"  Cluster,  and 
be  characteristic  of  an  even  younger  age. 

In  this  connection  we  should  note  that  the  unseen  com- 
panion  in  eclipsing  binary  BM  (=8'B)  Orionis,  hitherto  a 
very  puzzling  object  and  "black  hole"  candidate,  has  re- 
cently been  explained  by  Popper  (1975)  on  the  basis  of  his 
discovery  of  its  spectrum.   It  is  apparently  a  1 . 8M 
star  in  a  state  of  pre-main  sequence  differential  rotation. 


2  magnitudes  above  the  main  sequence.   The  currently  ac- 
cepted  contraction  time  for  such  an  object  is  ~1.8  •  10 
years . 

6.   Radio  observations  of  the  molecular  cloud  Orion  A 
have  recently  produced  two  surprises.   First,  the  mass  of 
material  (mostly  in  the  form  of  H^  and  not  observed  di- 
rectly but  deduced  from  the  observed  abundance  of  CO  and 

5 
HCN)  is  enormously  large:   about  10  Mg  (Liszt  et  al  ,  1974), 

or  several  hundred  times  that  of  the  stellar  mass  in  the 
same  volume.   If  the  nebula  really  is  expanding,  it  is  puz- 
zling why  so  large  a  fraction  of  the  material  failed  to 
form  into  stars . 

Second,  the  widths  of  the  CO  and  HCN  lines  are  much 
broader  than  the  thermal  Doppler  widths,  implying  variable 
mass  motion  along  the  line  of  sight.   Fitting  a  non-LTE 
model  to  the  data,  and  assuming  the  motions  to  be  large 
scale,  Gerola  and  Sofia  (1975)  are  led  to  conclude  a  linear 
contraction  amounting  to  12  km/sec  at  a  distance  of  1/2° 
from  the  center.   This  is  in  agreement  with  the  large 
gaseous  mass  of  the  cloud  and  apparent  range  of  stellar 
ages  but  not  with  the  idea  that  the  cluster  is  expanding. 
One  may  indeed  interpret  the  line  profiles  by  an  equal 
expansion  velocity,  but  then  one  finds  a  shell  midway  in 
the  cloud  where  the  gravitational  potential  energy  exceeds 
the  kinetic  energy. 


8 


On  the  other  hand,  Zuckerman  and  Evans  (197A)  have 
criticized  the  uniform  contraction  model,  proposing  instead 
that  the  observed  velocity  spread  is  due  to  turbulent  eddies 
interspersed  all  along  the  line  of  sight  through  the  cloud. 
If  this  is  so,  then  the  cloud  model  of  Gerola  and  Sofia  is 
invalid,  and  some  entirely  different  model,  as  yet  not 
worked  out,  would  be  required. 

The  Value  of  a  New  Astrometric  Study 
The  questions  fundamentally  at  issue  are  these: 

(1)  Is  the  Orion  Nebula  cluster  expanding,  as  is  the  Orion 
Association  as  a  whole,  or  is  it  still  undergoing  its  ini- 
tial contraction,  with  star  formation  still  taking  place? 

(2)  Are  the  motions  in  the  nebula  smoothly  varying,  or 
dominated  by  turbulence?   The  answers  to  these  questions 
are  important  to  our  understanding  of  the  process  of  star 
formation  in  clusters. 

A  new  study  of  the  kinematics  of  the  stars  In  the 
Orion  Nebula  would  form  a  valuable  basis  for  answering  these 
questions.   Thus  far,  the  astrometric  results  have  been 
discordant  and  inconclusive.   It  is  the  purpose  of  this 
investigation  to  derive  proper  motions  of  sufficiently  high 
accuracy  to  make  a  more  conclusive  kinematic  study 
possible . 


Scope  of  Present  Investigation 
Four  principal  desiderata  determined  the  scope  of  the 
present  study.   (1)  The  proper  motions  should  be  referred  to 
reference  stars  exclusively  and  In  no  way  depend  on  assump- 
tions about  the  optical  properties  of  the  instruments  with 
which  they  are  obtained.   We  call  motions  of  this  type 
"instrumentally  absolute."   In  particular,  they  are  to  be 
completely  independent  of  any  assumption  concerning  the 
variation  of  the  telescope's  plate  scale  with  time.   The 
proper  motions  are,  of  course,  relative  in  the  astrometric 
sense  in  that  they  are  referred  to  the  positions  and  mo- 
tions of  a  system  constituted  by  reference  stars  rather 
than  referred  to  the  equator  and  equinox  directly.   Within 
the  accuracy  of  the  determination,  the  positions  and  proper 
motions  derived  here  may  differ  from  the  best  derivable 
positions  and  motions  systematically  by  an  additive  con- 
stant, but  not  by  a  scale  factor.   (2)  The  coverage  should 
include  enough  reference  stars  to  determine  the  plate  scale 
as  accurately  as  possible.   (3)  As  many  stars  of  the 
Trapezium  cluster  as  practical  should  be  included,  while 
equally  faint  stars  farther  out  may  be  omitted.   (4)  The 
formulation  of  the  adjustment  equations  should  be  statisti- 
cally fully  rigorous;  that  is,  the  variance  of  the  residuals 
of  all  the  (weighted)  observations,  taken  together,  is  to  be 
minimized.   This  means,  among  other  things,  that  the  mathe- 
matics used  in  the  adjustment  must  take  into  account  the 
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fact  that  a  given  star  at  a  given  epoch  has  one  position 
regardless  of  from  whatever  plate  this  position  is  derived. 
Taking  due  account  of  this  latter  restraint  is  equivalent 
to  using  the  "overlap"  method  of  Eichhorn  (1963).   It  is 
the  imposition  of  the  overlap  condition  that  produces  not 
only  higher  internal  accuracy  than  otherwise  obtainable, 
but  also  allows  the  reference  stars  to  determine  more 
accurate  estimates  for  the  plate  scales  and  hence  the 
expansion.   This  is  so  because  there  are  not  enough  refer- 
ence stars  in  the  region,  and  therefore,  the  accuracy  of 
the  plate  scale  which  one  would  obtain  if  each  plate  were 
reduced  separately  from  all  others  is  insufficient  for  our 
purposes.   By  imposing  the  overlap  condition,  however,  we 
are  able  to  make  full  use  of  plates  on  which  as  few  as  no 
reference  stars  appear. 

A  fully  rigorous  solution,  in  the  strict  sense  defined 
above,  would  require  that  every  single  measurement  be  pre- 
served separately  throughout  the  reduction,  i.e.,  no  aver- 
aging of  individual  settings  of  the  measuring  engine,  or  of 
separate  grating  images,  could  be  made.   Since  strict  rigor 
is  laborious  to  enforce  and  really  unnecessary  when  the 
error  distribution  of  such  measurements  is  uncorrelated 
with  any  of  the  parameters  being  sought,  we  have  relaxed 
the  algorithm  to  allow  the  usual  practices  of  averaging 
direct  and  reverse  measures,  pairs  of  grating  images,  etc., 
in  the  interest  of  expediency. 


SECTION  II 
REFERENCE  MATERIAL 


By  reference  material  we  mean  the  positions  and  proper 
motions  of  the  set  of  reference  stars  with  respect  to  which 
the  positions  and  motions  of  all  the  program  stars  are  com- 
puted. Since  the  role  of  the  reference  material  is  espe- 
cially crucial  in  the  present  study,  it  deserves  discussion 
in  some  detail  here. 

As  usable  reference  material  we  may  only  include  posi- 
tions whose  errors  are  uncorrelated •   In  practice  such  posi- 
tions are  obtained  from  meridian  observations,  or  photo- 
graphic catalogues  sufficiently  global  that  our  selection 
of  reference  material  is  confined  to  a  relatively  small 
region.   Thus,  heliometer  or  micrometer  observations  are 
usable,  even  when  confined  to  a  small  region  of  the  sky;  so 
are  photographic  zone  catalogues,  if  the  zone  is  suffi- 
ciently wide,  even  though  such  positions  are  "secondary." 
On  the  other  hand,  V7e  consider  a  catalogue  based  on  photo- 
graphs covering  only  the  program  region,  e.g.  the  Zo-Se 
Catalogue,  unusable  as  reference  material  because  the 
errors  are  appreciably  correlated  through  the  plate  con- 
stant variances  of  the  photographs. 
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The  region  from  which  the  reference  material  is  taken 
is  approximately  bounded  by  5  26   <  a  <  5  42  ,  -6h°    <  6  < 
-4'5°  (orientation  1950).   The  actual  limits  were  dictated 
by  the  bounds  of  those  Astrographic  Catalogue  plates  which 
happen  to  include  the  program  region;  these  extend  beyond 
any  of  the  other  plate  material. 

Source  Catalogues.   We  have  taken  our  reference  mate- 
rial from  the  following  sources: 

1.  FK4  and  N30.   The  position  of  one  of  our  reference 

stars,  I  Ori,  is  listed  in  the  FK4 .   The  portions  of  three 

1     2 
additional  stars  are  found  in  the  N30:  of  6  C,  6  A,  and  of 

BD-6°  1255.   The  N30  positions  and  proper  motions  were 
transformed  to  the  FK4  system  according  to  the  tables  of 
Brosche  et  al .  (1964).   For  i  Ori,  the  mean  of  the  trans- 
formed position  and  that  of  the  FK4  has  been  used. 

2.  Boss'  General  Catalogue  (GC) ;   For  six  of  the 
reference  stars,  positions  are  found  in  only  the  GC.   The 
positions  and  motions  were  taken  directly  from  the  SAOC  and 
are  therefore  on  the  FK4  system. 

3.  Yale  Zone,  Volumes  16,  17:   The  bulk  of  the  refer- 
ence stars  is  listed  in  the  Yale  Zone  catalogues.   Their 
positions  and  motions  on  the  FK4  system  are  taken  from  the 
SAOC.   When  the  Yale  position  is  not  listed  in  the  SAOC,  it 
was  taken  from  vol.  1^  (Barney,  1945),  or  vol.  jL_7  (Barney, 
1945a),  and  transformed  to  the  FK4  system  by  adding  the 
corrections  (Yale  to  GC)  given  in  the  respective  introduction 
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and  the  correction  (GC  to  FK4)  from  Brosche,  et  al   (1964). 
Volumes  j^  and  JJ_  cover  the  zone  -10°  to  -2°,  joining  at 
-6°  with  vol.  1^  being  the  northern.   These  two  halves  had 
been  derived  from  one  single  belt  of  Yale  plates  and  re- 
duced as  one  zone;  the  proper  motions,  however,  were  ob- 
tained from  comparison  with  two  different  Strasbourg  AGKl 
zones  which  join  at  -6°.   A  comparison  of  the  proper  motions 
for  the  same  star  as  it  appears  in  both  volumes  for  24 
stars  in  common,  reveals  a  systematic  difference  between 
the  two  in  the  sense  y,  (vol.  Ij6)  =  p   (vol.  22)  -  V012/yr. 
We  have  therefore  added  the  correction  -'.'012/yr  to  all  the 
IJ^'s  from  vol.  j^.   If  this  is  not  done,  the  systematic 
difference  is  propagated  through  the  plate  constants  and  a 
marked  spurious  contraction  of  the  whole  region  in  6 
results . 

The  epoch  of  this  Yale  zone  is  1934.   The  AGKl  posi- 
tions (epoch  1895)  of  these  stars,  used  in  the  formation  of 
the  proper  motions  listed  in  these  Yale  volumes,  have  such 
large  errors  (a_V4) ,  that  the  resultant  error  in  the  proper 
motions  is  ,'.'01/yr,  corresponding  to  a  positional  error  of 
^'.'45  at  the  epochs  of  both  our  earliest  and  latest  plates 
(1900,  1974).   Accordingly,  it  is  very  desirable  to  improve 
the  proper  motions  of  these  reference  stars  with  additional 
reference  positions  of  epoch  well  before  or  after  1934. 

4.   Meyermann's  Heliometric  Positions;   B.  Meyermann 
(1903)  measured  the  positions  of  47  bright  stars  in  the 
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vicinity  of  the  Trapezium  in  1903,  using  the  Gottingen  heli- 
ometer.   This  series  is  especially  valuable  because  its 
early  epoch  and  the  relatively  high  accuracy  (_'.'25  in  each 
coordinate)  improves  the  accuracy  of  the  proper  motions 
given  in  the  Yale  Catalogue  from  VOl/yr  to  V007/yr. 
Thirty-seven  of  the  Meyermann  stars  are  included  in  the 
present  program;  of  these  the  data  for  31  are  in  the  Yale 
Catalogue;  for  five,  in  the  GC  only.   (The  remaining  star 
can  thus  be  used  as  a  reference  star  only  for  plates  taken 
near  the  epoch  1903.) 

Some  30  years  later,  in  the  course  of  a  photographic 
study  of  the  same  region,  Meyermann  discovered  an  error  in 
the  scale  of  the  heliometer;  accordingly,  in  a  subsequent 
paper  (Meyermann,  1938),  he  gave  revised  values  for  the 
1903  positions,  corrected  for  the  scale  error.   At  the  same 
time,  he  also  transformed  the  positions  to  the  system  of 
the  FK3  by  translation  and  rotation,  using  parameters  de- 
rived by  a  least-squares  fit  of  the  heliometer  positions  to 
those  in  the  PGC  and  the  "Courvois ler"   Catalogues.   The 
data  for  nine  stars  in  all  were  available  for  this  trans- 
formation.  The  relations  between  the  systems  of  the  PGC 
and  Courvoisier  catalogue  to  the  system  of  the  FK3  were 
known  and  used  by  Meyermann. 


Meyermann  does  not  identify  which  of  several 
Courvoisier  Catalogues  is  meant,  but  in  any  case  he  quotes 
the  actual  positions. 
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In  making  use  of  these  positions,  we  have  first  pre- 
cessed  the  coordinates  from  orientation  1900  to  1950  using 
Newcomb's  constants,  and  transformed  the  PGC  positions 
quoted  by  Meyermann  to  the  system  of  the  GC  (and  eventually 
to  the  FKA  system)  for  the  same  stars  by  a  least-squares 
adjustment  and  applied  that  transformation  to  all  the 
Meyermann  stars.   This  amounts  to  a  correction  Aa  =  -fo40, 
arising  from  a  systematic  error  in  the  PGC  in  the  amount  of 
Aa  =  -?060. 

5.   Southern  Reference  System  (SRS);   This  is  the 
reference  catalogue  currently  being  compiled  at  the  U.S. 
Naval  Observatory  from  relative  meridian  circle  observa- 
tions and  intended  as  a  counterpart  to  the  AGK3R  in  the 
southern  hemisphere.   It  contains  the  data  for  six  stars 
in  the  present  program,  all  of  these  included  also  in  the 
Yale  Catalogue.   The  average  epoch  is  1960,  thus  affording 
a  valuable  improvement  in  the  accuracy  of  their  proper 
mo  t  ions . 

Preliminary  positions  for  these  stars  were  made  avail- 
able prior  to  publication  of  the  completed  catalogue 
through  the  courtesy  of  Dr.  J.  Schombert  of  the  USNO. 

Combination  of  reference  sources.   In  the  case  of 
stars  Included  in  more  than  one  of  the  reference  sources, 
a  weighted  mean  position,  proper  motion,  and  epoch  has  been 
computed  according  to  the  rigorous  expressions  (eg. 
Eichhorn,  1974   p.  113).   An  exception  is  made  in  the  case 
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of  the  three  FK4  and  N30  stars,  these  positions  not  being 
averaged  with  those  from  other  sources.   The  standard  devi- 
ations characteristic  of  each  source  are  listed  in 
Table  1. 

Double  Stars.   The  Aitken  Double  Star  Catalogue  (ADS) 
was  consulted  for  any  entries  lying  within  the  range  of  our 
program.   None  of  our  program  stars  shows  any  clear  sign  of 
orbital  motion  from  the  observations  listed  in  the  ADS,  so 
no  orbital  corrections  have  been  applied.   Two  uses  were 
made  of  the  ADS  material: 

1.  Doubles  closer  than  3"  but  wider  than  1"  were  re- 
jected from  the  reference  list  since  they  might  appear 
variously  resolved  on  some  plates  and  unresolved  on  others. 
Some  wider  doubles  (up  to  5")  were  also  rejected  on  the 
shorter-focus  astrographic  plates  for  the  same  reason. 
Doubles  closer  than  1"  have  all  been  retained  since  they 
are  presumably  always  unresolved.   These  criteria  lead  to 
the  rejection  of  the  following  reference  stars:   i  Ori  (on 
AC  plates),  BD  -4.1171,  -4.1172,  -4.1185,  -4.1186  (on  AC 
plates),  and  -6.1255  (on  AC  plates). 

2.  Astrometric  data  for  doubles  were  utilized  in  the 
following  special  cases:   (a)  Wide  pairs  which  occur  as  two 
distinct  reference  stars,  all  of  which  happen  to  show  zero 
relative  motion,  were  assigned  equal  reference  proper  motion 
to  both  members,  formed  by  the  mean  of  the  two  catalogue 
values.   This  amounted  in  some  cases  to  a  substantial 


TABLE  1 
SUMMARY  OF  REFERENCE  CATALOGUE  MATERIAL 


19 


Catalogue 


a  ,  s 


FK4 

V05 

N30 

.05 

Yale  (v.  16  or  17) 

.14 

GC 

.15 

Meyerraann  (H) 

.25 

S.R.S. 

.05 

Yale  (v.  16  &  17) 

.10 

Yale  +  GC 

.10 

Yale  +  H 

.12 

Yale  +  S.R.S. 

.05 

V004/yr 

1912 

.004/yr 

1912 

.010/yr 

1934 

•OlO/yr 

1900 

.007/yr 

1903 

.007/yr 

1965 

.007/yr 

1934 

.007/yr 

1915 

.007/yr 

1924 

.003/yr 

1960 
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departure  from  the  original  catalogue  value,  and  always 
yielded  a  better  fit  to  the  plate  measurements.   This  pro- 
cedure does  not,  of  course,  constrain  the  final  proper 

motions  of  the  components  to  be  equal.   The  following  pairs 

1        2     2        2 
were  treated  this  way:   e  C  and  6  A;  6  A  and  9  B;  BD 

-6°1233  and  1234.   (b)  For  S'^'E  and  S-'t,  the  fainter  com- 
panions to  the  four  principal  Trapezium  stars,  the  micro- 
metric  observations  prior  to  1920  are  used  in  the  deter- 
mination of  their  proper  motions,  since  these  two  stars 
occur  on  no  early  epoch  plates. 

Parallaxes .   Only  two  of  our  program  stars  have  known 
and  measurable  parallaxes:   i  Ori  (tt  =  V02)  and  45  Orl  (tt  = 
'.'02)  .   These  values  are  too  small  to  influence  the  overlap 
solution,  but  are  probably  large  enough  to  establish  the 
stars  as  foreground  objects.   Correction  for  the  effects  of 
these  parallaxes  have  been'applied  to  their  reference 
posit  ions . 

Magnitudes  and  colors.   In  general,   we  must  antici- 
pate that  the  plate  reduction  models  may  include  terms  in- 
volving magnitude  and  color.   Therefore,  we  need  approxi- 
mate values  for  them.   For  this  purpose,  they  need  to  be 
known  with  a  relative  accuracy  of  about  0™1.   As  a 
criterion  for  membership  in  the  cluster,  however,  we  re- 
quire more  accurate  values.   As  data  we  have  (1)  used  the 
photoelectric  photometry  of  Johnson  (1957)  and  Sharpless 

(1962)  when  available,  and  (2)  computed  m„  and  m_  for  the 

V        B 
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remaining  stars  from  measurements  on  the  iris  photometer 
of  two  pairs  of  blue  and  yellow  plates,  each  pair  taken 
simultaneously  especially  for  the  purpose,  and  using  the 
Johnson  and  Sharpless  stars  as  standards.   The  details  and 
final  results  of  this  photometry  will  be  published  else- 
where; the  colors  given  here  should  be  considered 
preliminary . 


SECTION  III 
THE  PLATE  MATERIAL 


General  Survey 
In  this  sub-section  the  whole  of  the  plate  material 
is  summarized.   The  next  two  sub-sections  deal  with  the  USF 
plates  particularly  and  with  the  plate  measurements  carried 
out  in  the  present  program. 

The  plate  material  (that  is,  the  x,y  coordinates  of 
the  images  of  the  program  stars  as  measured  on  the  photo- 
graphic plates)  used  in  this  study  is  very  heterogeneous. 
This  is  so  partly  by  necessity  and  partly  by  design.   There 
exists  no  one  series  of  plates,  taken  at  one  telescope, 
having  at  once  all  four  of  these  properties:   (a)  epoch 
difference  of  at  least  30  years,  (b)  plate  scale  of  at 
least  (say)  20y/",  (c)  inclusion  of  an  adequate  number  of 
reference  stars,  and  (d)  measurable  images  of  the  inner 
cluster  stars  (i.e.  stars  within  5'  of  the  Trapezium). 
Thus,  the  Yerk.es  plates,  for  example,  bear  images  of  only 
seven  reference  stars;  the  Allegheny  plates  are  obscured 
at  the  center  by  over-exposed  nebulosity,  while  on  the 
McCormick  plates  the  central  region  lacks  all  but  the 
brightest  stars  due  to  the  use  of  a  sector  wheel.   The 
U.S.F.  and  Sproul  plates  avoid  all  these  difficulties,  but 
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have  no  early  epoch  counterparts.   As  for  the  possibility 
of  repeating  the  AC  type  "as trographic  camera"  photography 
of  this  region  with  modern  plates  from  the  same  telescopes, 
this  is  not  feasible  in  the  present  instance  because  the 
San  Fernando  plates  covering  the  region  had  an  exception- 
ally bright  limiting  magnitude  (^11  ),  while  the  early 
epoch  "Zo-Se"plates  (of  excellent  quality)  cannot  be  re- 
peated because  the  telescope  with  which  they  were  taken  has 
disappeared.   In  any  case,  the  focal  length  of  the  astro- 
graphic  cameras  is  woefully  short — about  one-third  that  of 
a  typical  long  -focus  "parallax"  instrument  such  as  the 
McCormick.  refractor. 

It  is  regrettable  that  an  objective  grating  was  not 
employed  in  any  of  the  early  epoch  photography. 

The  plate  material  is  summarized  in  Table  2.   Explana- 
tory notes  and  references  are  provided  in  the  caption.   A 
list  of  the  plate  numbers  of  the  individual  plates  used,  in 
the  notation  of  their  respective  observatories,  is  given  in 
the  Appendix  (Table  A-1).   The  information  concerning  the 
circumstances  of  the  USF  photography,  since  it  has  not 
appeared  elsewhere,  is  also  listed  in  the  Appendix 
(Table  A-2) . 

The  USF  Photography 


Telescope 

A  description  of  the  astrometric  reflector  of  the 
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TABLE  2 
SUMMARY  OF  PLATE  MATERIAL 


Plate 

No  .  o 

f 

Scale 

Epoch 

No.  of 

No.  of 

Ref  . 

d 

Q 

f 

Obs^ 

(y/") 

Plates 

Starsc 

Stars 

Msmt 

Source 

S  F 

17 

1893,1906 

3 

130 

85 

M 

1 

Z-S 

17 

1904-1916 

5 

M 

2 

Y 

93 

1905-1912 

7 

92 

7 

M 

3 

A 

68 

1921 

6 

70 

15 

M,A 

4 

McC 

48 

1924-1927 

22 

47 

13 

A 

5 

McC 

48 

1968 

22 

47 

13 

A 

5 

A 

68 

1970 

5 

70 

15 

A,M 

4 

Sp 

53 

1974 

2 

90 

14 

A 

4 

USF 

50 

1974 

20 

140 

35 

A,F 

4 

^Name  of  observatory:  S  F  =  San  Fernando,  Spain;  Zo-Se, 

also  "Zi-Ka-Wei,"  Shanghai,  China;  Y  =  Yerkes ;  McC  =  McCormick; 

A  =  Allegheny;  Sp  =  Sproul;  USF  =  Observatory  of  University  of 
South  Florida,  Tampa. 

Epoch:   approximate  mean  value,  or  range. 

'^No.  of  stars:   number  of  program  stars  appearing  on  at 
least  2  plates. 

No.  of  ref.  stars:   number  of  reference  stars  with  mea- 
surable images  on  at  least  2  plates.   (The  components  of  6 
and  9^  are  counted  here  as  only  two  ref.  stars.) 

^Msmt:   Type  of  measuring  engine--M  manual;  A  automatic 
(impersonal);  F  fringe- count ing  (screw- independent ) . 

Source:   source  of  measurements  as  follows: 

(1)  Astrographic  Catalogue,  San  Fernando  Section,  4  (-5° 
zone),  p.  50  (plate  1496);  and  5  (-6°  zone),  p.  50  (plates 
1490,  3890). 

(2)  Chevalier  (1933).   The  listing  for  plate  "2"  con- 
tains many  star  misidentif ications  discovered  and,  where  pos- 
sible, rectified  by  this  author.   There  are  also  several  typo- 
graphical errors  in  the  measurements.   This  photography  has 

no  relation  to  the  Zo-Se  equatorial  zone,  published  shortly 
before  in  the  same  series  of  publications.   In  general,  all 
of  the  photography  and  measurements  carried  out  under  the 
supervision  of  Chevalier  at  Zo-Se  were  meticulously  and  in- 
geniously done  so  as  to  yield  very  high  accuracy. 

(3)  Strand  (1958,  1972)  for  description;  the  measure- 
ments on  punch  cards  provided  by  Dr.  Strand. 
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Table  2--contlnued 

(4)  Present  investigation.   Measurements  made  in  1973 
and  1974,  described  below. 

(5)  Cannell  (1970)  for  a  description;  the  measurements 
on  punched  cards  provided  by  Mr.  Cannell.   Both  the  early 
and  late  epoch  plates  were  measured  by  him  in  1970. 


University  of  South  Florida  (USF)  Observatory  has  not  yet 
appeared  in  the  literature.   The  telescope  is  a  26-inch 
Tinsley,  f/15  "Schmidt-Cas segrain"  with  optics  designed  by 
J.  G.  Baker.   These  are  a  full-aperture  corrector  plate 
(not  a  "reflector- corrector" ), an  f/5   nearly  spherical 
CERVIT  primary,  and  an  approximately  spherical  secondary  in 
a  Cassegrain  configuration.   This  system  has  also  been 
termed  a  "Baker-Schmidt."   The  unobstructed  field  is  1?2. 
Stellar  images  are  sharp  and  round  up  to  the  edge  of  the 
field,  exhibiting  no  trace  of  coma. 

As  is  the  case  with  systems  of  the  Cassegrain  type, 
image  quality  is  critically  sensitive  to  the  collimation 
of  the  elements.   With  the  USF  telescope,  proper  collima- 
tion is  achieved  only  with  difficulty;  moreover,  the  mirror 
supports  are  such  that  collimation  cannot  be  maintained  for 
more  than  a  few  weeks  at  a  time.   Consequently,  the  condi- 
tions under  which  the  plates  used  for  this  investigation 
were  obtained  cannot  with  certainty  be  duplicated  at  a 
later  time,  and  such  instrumental  plate  constants  which 
might  otherwise  have  remained  constant  over  time  cannot  be 
assumed  to  have  done  so  here. 
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The  USF  Observatory  Is  located  within  an  urban  area 
and  suffers  from  a  sky  brightness  that  sets  the  limiting 
magnitude  at  about  14   on  yellow  plates. 

Coverage 
The  format  of  the  USF  plates  is  5  x  7  inches.   Cen- 
tered at  the  Trapezium,  this  format  encompasses  only  about 
one  dozen  reference  stars.   In  order  to  utilize  more  ref- 
erence stars  we  must  extend  the  field  by  overlapping  adja- 
cent plates.   The  overlap  pattern  chosen,  shown  in  Figure  2, 
was  dictated  by  both  the  distribution  of  the  reference 
stars — which  is  far  from  uniform--and  limitations  of  eco- 
nomy.  It  is  seen  that  there  are  very  few  reference  stars 
east  and  especially  west  of  the  cluster--a  circular  area 
one  full  degree  in  diameter  centered  about  1/2°  west  of  the 
Trapezium  contains  but  two  reference  stars.   The  distribu- 
tion of  plates  is  center-to-corner,  with  the  Trapezium  ap- 
pearing on  nearly  all  the  plates;  the  region  covered  is 
roughly  co-extensive  with  the  "Sword"  of  Orion,  encompass- 
ing a  total  of  39  reference  stars.   Sixteen  of  these  plates 
were  used  in  this  study. 


Exposures 
Photography  of  the  Orion  Trapezium  cluster  presents 
several  challenges.   First,  of  course,  is  the  nebulosity 
itself.   On  the  plate,  the  background  light  of  the  nebula 
blots  out  the  fainter  stars  imbedded  in  it,  while 
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condensations  of  the  nebula  adjacent  to  stellar  images  can 
shift  the  apparent  photo-centers  of  those  images.   Second, 
the  magnitude  range  of  the  program  stars  is  at  least  seven 
magnitudes.   Within  the  vicinity  of  the  Trapezium  espe- 
cially, bright  and  faint  stars  are  closely  intermingled, 
precluding  the  use  of  a  sector  wheel.   An  objective  grating 
can  only  be  used  with  great  care,  lest  the  multitude  of 
images  of  the  different  stars  impinge  upon  one  another. 

We  have  sought  to  circumvent  these  difficulties  with 

the  following: 

1.   Choice  of  filter  and  emulsion.   All  exposures  were 
made  on  103a-G  emulsion,  taken  through  a  Corning  3-69  fil- 

o 

ter.   The  resultant  combination  isolates  a  600  A-wide 
region  free  of  major  nebular  emission  lines.   The  suit- 
ability of  this  combination  is  revealed  in  Figure  1,  in 
which  the  nebulosity  is  clearly  suppressed. 

The  filter  can  be  expected  to  introduce  some  distor- 
tion of  the  field.   In  order  to  counteract  this,  we  posi- 
tion the  filter  immediately  in  front  of  the  plate,  accord- 
ing to  standard  practice,  the  spacing  being  3ram.   In  addi- 
tion, we  have  taken  the  unusual  precaution  of  calibrating 
the  filter,  as  described  below  under  "Measurement  Process- 
ing."  To  facilitate  in  applying  this  calibration,  a  ref- 
erence mark  is  scribed  on  the  filter  so  as  to  be  sil- 
houeted  on  plates  exposed  through  it. 
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2.  Objective  grating.   A  coarse  26-bar  objective 
grating  was  constructed  and  employed  in  some  of  the  expo- 
sures.  The  direction  of  the  bars  was  rotated  until  the 
diffraction  patterns  of  the  Trapezium  stars  all  avoided  one 
another — this  proving  a  most  delicate  adjustment.   The 
grating  constants  were  found  by  experiment  to  be  4.4,  4.7, 
5.0,  6.3,  respectively  for  the  first-through  4th-order 
images.   (In  retrospect,  it  seems  probable  that  a  more 
evenly  graded  set  of  constants  would  have  been  more  useful.) 

3.  Multiple  exposures.   As  an  alternative  to  the  use 
of  the  grating,  some  plates  were  taken  with  both  a  long  and 
a  short  exposure.   An  x,y  translation  in  the  focal  plane  of 
less  than  5  mm  was  found  such  that  the  images  of  one  expo- 
sure would  not  fall  on  any  of  those  of  the  other.   Exposure 
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times  or  10     and  1     were  usually  made,  separated  in  time 

by  only  a  few  seconds  while  the  plate  holder  was  advanced 

in  the  camera  by  the  pre-set  distance.   Some  of  the  plates 

taken  with  the  grating  were  also  given  two  exposures. 

Photographs  were  taken  on  either  side  of  the  pier, 
usually  at  less  than  30  min.  from  the  meridian  but  occa- 
sionally more.   Exposure  times  were  usually  9  minutes, 
reaching  IZh    .   Particulars  are  recorded  in  Table  A-2. 

All  plates  were  kept  horizontal  during  developing  and 
drying  to  reduce  emulsion  shifts.   For  the  same  reason, 
washing  time  was  reduced  to  a  bare  minimum  even  at  the  ex- 
pense of  incurring  cosmetic  blemishes,  while  drying  time 


29 

was  reduced  by  using  a  "drying  box"  with  circulating  warm 
air. 

Measuring  the  Plates 

General  Consideration 

As  noted  above,  measurements  for  some  of  the  plates 
used  in  the  present  study  were  already  available.   The 
other  plates  were  measured  in  the  course  of  this  investiga- 
tion.  The  latter  are  20  USF  plates,  12  Allegheny  plates, 
and  two  Sproul  plates. 

Ideally,  one  of  the  automatic  measuring  engines,  such 
as  the  SAMM  (Strand  Automatic  Measuring  Machine)  of  the 
U.S.  Naval  Observatory,  would  have  been  used  for  the  mea- 
surement of  all  plates  because  of  both  its  impersonal  cen- 
tering action  and  much  greater  speed.   In  practice,  we 
found  it  is  not  always  possible  to  use  such  a  machine,  for 
the  following  reasons.   (The  limitations  referred  to  here 
apply  specifically  to  the  SAMM,  but  similar  limitations 
apply  as  well  to  the  other  automatic  machines.)   The  effec- 
tive scanning  area  perceived  by  the  SAMM  centering  mechan- 
ism is  a  circle  of  250m  diameter.   On  the  SAMM  it  is  impos- 
sible to  properly  measure  an  image  which  extends  beyond 
this  circle,  or  which  is  not  isolated  within  the  circle 
which  happens,  for  instance,  when  a  part  of  another  image 
also  falls  within  the  circle;  or  when  the  background  varies 
measurably  across  the  circle's  area  or  is  darker  than 
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density  .5.   These  limitations  are  often  exceeded  in  the 
Orion  nebula  photography,  due  to  the  crowding  together  of 
grating  images,  the  presence  of  the  nebulosity  itself,  or 
(in  the  case  of  some  of  the  USF  photography),  the  presence 
of  sky  fogging.   The  large  range  of  magnitudes  encountered 
in  our  program  also  has  the  result  that,  on  plates  taken 
without  a  grating,  the  brighter  stars  (these  include  nearly 
all  the  reference  stars,  which  must  be  measured)  produce 
images  with  diameters  larger  than  the  critical  250^.   This 
phenomenon  is  especially  at  its  worst  on  the  early 
Allegheny  plates  which  were  taken  without  the  filter  re- 
quired to  compensate  for  the  difference  in  spectral  regions 
for  which  the  chromatic  aberration  of  the  objective  is  cor- 
rected and  that  to  which  the  plates  were  sensitive. 

The  plates  were  inspected  visually  to  judge  the  degree 
to  which  these  effects  were  present  and  each  plate  was 
accordingly  assigned  to  be  measured  either  on  the  SAMM  or 
on  a  manual  machine.   Some  plates  were  measured  in  both 
modes-~automat ic  for  the  smaller  images  and  manual  for  all 
the  images.   A  brief  account  of  the  respective  measuring 
procedures  follows. 


Measuring  Procedures 

The  automatic  machine.   The  SAMM  of  the  U.S.  Naval 
Observatory  was  made  available  through  the  courtesy  of 
Dr.  P.  M.  Routly  in  1973  and  in  1974  and  used  for  the 
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measurement  of  all  these  plates  which  could  be  measured  on 
it.   This  machine  has  been  described  elsewhere  (Strand, 
1971).   Measurement  sets  were  made  in  "direct"  and  "re- 
verse" orientations  and  combined  by  a  least-squares  adjust- 
ment (hereafter  called  a  "rotation")  solving  for  transla- 
tion, rotation,  and  a  magnitude  term,  just  as  is  conven- 
tionally done  for  manually  made  measurements.   The  resid- 
uals (direct  minus  final)  from  these  adjustments  have  dis- 
persions of  about  1.5  microns. 

The  calibration  corrections  given  by  the  USNO  for  this 
machine  were  applied  to  the  measurements  made  in  1973.   It 
was  found  that  the  measurements  made  in  1974,  after  a  series 
of  adjustments  had  been  made  to  the  SAMM,  were  better 
fitted  direct  to  reverse  without  the  corrections,  so  the 
corrections  were  not  applied  to  them.   At  most,  the  cor- 
rections amount  to  2  microns. 

The  "fringe-counting"  machine.   Also  at  the  USNO, 
there  is  a  conventional  Mann  two-screw  measuring  machine 
which  has  been  modified  so  that  the  stage  travel  is  mea- 
sured by  Moire  fringe-counters,  as  in  many  automatic  ma- 
chines, while  the  transport  and  centering  remain  unchanged 
(hence  this  machine  is  called  a  "screw-independent"  manual 
machine).   About  half  of  the  USF  plates  were  measured  on 
this  machine;  the  Allegheny  plates  are  too  large  to  fit  on 
its  stage.   The  temperature  at  the  machine  could  be  mea- 
sured but  not  held  constant;  therefore,  measuring  was 
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halted  when  the  temperature  varied  by  more  than  1°5C.   A 
calibration  of  the  machine  that  had  been  made  at  the  USNO 
shows  that  the  corrections  are  smaller  than  one  micron, 
therefore,  none  were  applied.   Four  settings  were  made  on 
each  image  in  both  0'^'  and  180°  orientations,  and  the  two 
orientations  combined  as  described  above.   After  the  fit, 
the  rms  error  of  the  effective  combined  measurements  is 
1.75p. 

Manual  measuring  machine.   A  conventional  two— screw 
Mann  engine  at  USF  was  used  to  measure  most  of  the  Allegheny- 
places  and  a  few  of  the  USF  plates.   The  temperature  of  the 
machine  remained  within  a  1°C  range  during  operation.   Four 
settings  per  image  were  made  in  each  of  four  orientations: 
0°,  90°,  180°,  and  270°,  thus  giving  two  pairs  of  direct 
and  reverse  measurement  sets.   Each  direct-reverse  pair  was 
"rotated"  as  above,  then  the  two  resulting  combined  sets 
also  rotated  one  into  the  other.   The  dispersion  of  the 
residuals  after  adjustment  confirms  the  value  of  measuring 
four  orientations:   for  each  initial  rotation  the  rms  error 
is  about  2h    microns,  while  for  the  final  rotation,  it  is 
1.5  microns.   Thus  the  two-fold  increase  of  information 
contained  in  four  as  opposed  to  two  orientations  is  com- 
pletely effective  in  reducing  the  mean  error  of  the  final 
combined  measured  coordinates. 

Automatic  and  manual  machines  used  jointly.   The 
plates  measured  on  both  the  SAMM  and  manual  machines  were 
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treated  as  follows.   At  the  SAMM,  all  measurable  Images 
were  measured,  in  direct  and  reverse  positions,  and  the 
measurement  sets  rotated  and  combined  as  for  the  other 
plates.   At  the  USF  machine  all  the  images  were  measured, 
and  the  measurements  were  likewise  rotated  as  above.   Then 
the  combined  "manual"  set  was  rotated  to  fit  the  corres- 
ponding "automatic"  set,  using  as  the  transformation  terms 
translation,  rotation,  and  scale.   The  rms  error  of  this 
adjustment  was  typically  2\i  .       Finally,  the  transformation 
parameters  were  applied  to  the  remaining  manual  measure- 
ments, producing  a  presumably  homogeneous  set  of  measure- 
ments for  all  the  images.   We  recognize  that  this  procedure 
of  combining  the  measurements  from  two  machines  is  not 
fully  rigorous.   This  is  so  because  the  adjustment  para- 
meters themselves  contain  some  error,  and  this  error,  which 
would  not  be  present  if  the  sets  were  not  combined,  infects 
the  manually-made  measurements.   The  rigorous  procedure 
would  be  to  segregate  the  two  sets,  and  during  the  plate 
reductions  allow  separate  translation,  rotation  and  scale 
terms  to  be  associated  with  each  set.   Since,  however,  at 
least  20  stars  common  to  both  the  SAMM  and  the  manual  mea- 
surements, distributed  uniformly  over  the  plate,  were 
available  in  all  cases  for  the  adjustment,  we  feel  no 
errors  greater  than  1  micron  have  been  introduced  by  hav- 
ing adopted  this  procedure. 
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Preliminary  Processing 

It  will  suffice  here  to  give  just  a  brief  outline  of 
the  subsequent  processing  of  the  measurements.   We  have 
taken  into  consideration  the  following: 

Filter  calibration.   It  should  be  anticipated  that  a 
filter  may  shift  the  positions  of  the  images  as  they  appear 
on  the  plate,  since  the  filter  may  well  depart  enough  from 
optical  flatness  to  introduce  a  measurable  and  non-uniform 
distortion  of  the  field.   Such  a  distortion  would  not  ap- 
pear as  a  degradation  in  image  quality,  nor  necessarily  be 
revealed  at  any  step  in  the  reductions.   We  control  the 
problem  by  "calibrating"  the  filters  used  at  the  USF  Ob- 
servatory.  This  is  done  by  placing  on  the  stage  of  the 
measuring  engine  a  grid  and  over  it  the  filter  in  exactly 
the  same  configuration  as  the  photographic  plate  and  filter 
occur  at  the  telescope,  then  measuring  the  coordinates  of 
the  grid  line  intersections  with  and  without  the  filter  in 
place.   The  shifts  in  the  coordinates  of  the  grid  then  pre- 
sumably duplicate  the  distortion  produced  by  the  filter  at 
the  telescope.   The  corrections  are  easily  applied  by  mea- 
suring the  image  of  the  filter  reference  mark  on  the  plates 
(see  "USF  Photography"  above)  at  the  same  time  the  star 
images  are  measured.   The  filter  "Y"  used  for  the  present 
USF  photography  proved  to  require  no  corrections  even  as 
large  as  1  micron. 
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No  calibration  has  been  made  of  the  filters  used  for 
the  other  plates  in  our  study,  i.e.  the  McCormlck,  later 
Allegheny  and  Sproul  plates.   It  is  recognized  that  system- 
atic errors  of  position  and  motion  may  remain  by  failing  to 
have  done  so,  especially  if  different  filters  were  used  for 
the  early  and  late  plates  at  the  same  telescope.   In  any 
case,  the  filter  used  for  the  early  McCormick  plates  was 
broken  shortly  after  the  plates  were  taken. 

Grating  images.   Where  higher-order  diffraction  images 
appear,  the  mean  was  taken  of  each  pair  of  corresponding 
images  and  subsequently  treated  as  a  single  meafjurement  for 
that  star,  with  its  magnitude  taken  as  the  magnitude  of  the 
star  plus  the  appropriate  grating  constant.   It  can  be 
shown  that  second-order  non-linearity  of  the  plate  scale 
over  a  one-degree  field  introduces  through  this  practice  an 
error  of  less  than  .Olp.   The  various  orders  were  not 
averaged  together. 

The  dispersion  of  the  differences  between  the  central 
image  and/or  means  of  the  higher-order  image  pairs  is  3.3y 
for  both  X  and  Y  for  the  USF  plates  and  4 . Ip  for  the  Sproul 
plates . 

Multiple  exposures.   For  plates  with  two  (or  three) 
exposures,  the  average  translation  (AX, AY)  between  the 
first  and  the  second  (and  third)  exposures  was  found.   A 
linear  dependence  of  (AX, AY)  on  magnitude,  as  might  arise 
from  a  difference  in  magnitude  term  between  exposures,  was 
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allowed  for.   The  values  (  X,  Y)  were  applied  to  all  the 
second  (and  third)  exposure  images,  whether  or  not  a  first 
exposure  image  was  present,  to  form  one  homogeneous  set  of 
measurements.   If  a  small  rotation  between  the  exposures 
also  is  present,  that  rotation  will  be  incorporated  into 
the  final  meaned  coordinates  with  no  detrimental  result; 
however,  a  systematic  rotation  of  the  set  of  single  expo- 
sures with  respect  to  the  set  of  means  of  two  exposures  will 
be  introduced.   This  procedure  is  rigorously  valid  if  only 
the  plate  was  moved  between  exposures.   If  the  tangential 
point  was  also  shifted  (as  in  re-centering  the  guide  star) 
an  error  is  introduced  amounting  to  l|i  over  a  1°  field  only 
when  the  shift  exceeds  5  mm. 

On  two  USF  plates,  the  shorter  exposures  appeared  to 
be  somewhat  elongated.   In  these  cases,  the  separate  expo- 
sures were  treated  as  though  they  were  separate  plates. 
Several  McCormick  plates  also  bore  two  exposures  each;  these 
were  all  treated  as  separate  plates. 

Refraction  and  aberration.   Since  the  apparent  posi- 
tions of  the  reference  stars  are  altered  by  refraction  and 
aberration  just  as  are  the  positions  of  the  other  program 
stars,  it  is  only  the  differential  values  of  these  effects, 
across  our  2°  region,  which  need  concern  us.   Moreover, 
since  our  plate  reduction  models  include  linear  and  quad- 
ratic terms  with  coefficients  to  be  determined  separately 
for  X  and  Y,  it  is  only  third-order  differential  refraction 
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and  aberration  which  remain  explicity  to  be  allowed 
for. 

The  largest  zenith  distance  at  which  any  of  the  plates 
was  taken  is  45°  .   At  this  value  the  maximum  second-order 
differential  refraction  is  '.'016;  the  third-order  differen- 
tial is  less  than  '.'001.   Third-order  differential  abbera- 
tion  is  also  less  than  '.'001.   Consequently,  corrections  for 
refraction  and  aberration  have  not  been  applied. 


SECTION  IV 
FORMULATION  OF  THE  PROBLEM 


The  method  used  for  the  computation  of  positions  and 
motions  in  this  study  falls  into  that  general  category 
termed  "plate  overlap"  methods.   In  presenting  our  formula- 
tion, we  will  develop  it  from  the  context  of  earlier  me- 
thods, beginning  with  the  simplest,  so  that  points  of  de- 
parture are  more  clearly  seen. 

Differential  Plate  Measurements 
The  most  direct  computation  of  proper  motion  is  simply 
the  comparison  of  the  coordinates  of  the  same  star  as  it 
appears  on  plates  taken  with  the  same  telescope  at  differ- 
ent epochs.   The  differences,  after  terms  in  arbitrary 
translation  rotation,  and  scale  have  been  removed  by  a 
least-squares  adjustment,  yield  the  motions  directly.   We 
may  write  this  formulation  as  follows.   Let  (x,y)  be  the 
coordinates  of  a  star  on  any  plate,  and  x  ,  y   be  initial 
coordinates  of  the  star,  say  the  coordinates  from  one  of 
the  plates  chosen  as  a  standard,  at  t  =  0.   For  convenience 

we  can  write  u  ,  li   for  the  components  of  u    parallel  to  the 
X   y 

axes  of  X  ,  y  ,  multiplied  by  the  focal  length.   Then  the 
oo 

coordinates  and  motions  for  a  star  on  the  m-th  plate  are 
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related  by  observation  equations  involving  a  function  f   of 

m 

parameters  (p  )  and  epochs  t   for  that  plate: 

f    (x,y,p)     =    (p    +p-x+p,y)       =    X   +(yt    +Ax)     e    x   +r    .  (1) 

ID  izjtnom  om 


For  brevity,  we  let  f   stand  for  the  vector  of  two  equa- 
tions,  and  n,  x  ,  Ax,  and  r  stand  respectively  for  the 
vectors  (yx,My),  (x  ,y  ),  etc.   The  parameters  (p  )  are 
found  by  a  least-squares  adjustment  for  each  of  the  M 
plates.   The  term  in  parentheses  on  the  right  side  repre- 
sents the  "residual"  for  each  star;  it  contains  a  random 
error  component  and  a  component  due  to  the  proper  motion. 
The  computation  of  the  y  from  the  M  residuals  can  be  done 
rigorously  by  writing 


|j  t   =  r 
m     m 


(2) 


and  solving  for  y  by  a  least-squares  adjustment.   It  is 
more  customary  to  combine  the  material  in  pairs  of  early 
and  late  plates  and  compute  the  mean  of  the  individual  prop- 
er motions,  weighted  according  to  the  individual  time 
interval.   It  is  to  be  noted  that  this  will  not  yield 
exactly  the  same  value  for  p  as  will  the  least-squares  re- 
duction, unless  the  ranges  of  early  and  late  epochs  sepa- 
rately are  small  compared  to  the  whole  range  of  epochs. 
The  great  advantage  of  this  method  (apart  from  its 
simplicity)  is  that  as  long  as  the  plate  centers  are 
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approximately  the  same  (a  condition  usually  met  in  prac- 
tice) ,  the  higher-order  terms  in  the  imaging  properties  of 
the  telescope  do  not  appear.   This  means  that  the  corres- 
ponding errors  in  those  terms  do  not  appear  in  the  errors 
of  the  proper  motions;  and,  since  such  plate  constant 
errors  increase  rapidly  toward  the  edge  of  the  plate 
(Eichhorn  and  Williams,  1963),  the  proper  motions  as  de- 
rived by  plate  differences  are  equally  accurate  out  to  the 
edge  of  the  field,  at  least  in  comparison  with  proper  mo- 
tions derived  through  classical  plate  reduction.   (Of 
course,  one  may  formulate  the  problem  this  way  even  when 
the  plates  have  been  obtained  with  different  telescopes, 
but  then  the  need  for  higher-order  terms  in  the  observation 
equations  vitiates  the  main  advantage  of  the  method.) 

The  disadvantage  is  that  the  terms  in  translation, 
rotation,  and  scale  absorb  any  net  constant  proper  motion, 
rotation  or  expansion  of  the  stellar  group  itself.   An  at- 
tempt is  sometimes  made  to  pin  down  these  terms  by  imposing 
external  constraints:   the  change  in  plate  scale,  for  exam- 
ple, can  be  monitored  by  a  field  of  known  proper  motions, 
photographed  at  the  epochs  of  the  earliest  and  latest 
plates.   Again,  the  rotation  might  be  monitored  by  trailing 
some  of  the  plates.   The  trouble  with  these  devices  is  that 
they  do  not  really  eliminate  the  problem.   As  long  as  any 
factors  which  might  give  rise  to  rotation  and  scale  terms, 
such  as  differences  in  refraction,  plate  tilt,  measuring 
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engine  parameters,  etc.,  remain  unknown,  adjustment  para- 
meters for  them  must  still  be  included  in  the  reduction. 

An  example  of  the  application  of  this  method  is 
Strand's  (1958)  study  of  the  Trapezium  cluster.   As  ex- 
pected, the  internal  dispersion  of  proper  motions  is  very 
small  (except  as  noted  above),  while  his  value  for  the 
expansion  has  been  subject  to  question. 

Reference  Proper  Motions 
The  disadvantages  mentioned  in  the  preceding  section 
may  be  eliminated  while  retaining  the  advantages  of  a  dif- 
ferential method,  by  making  use  of  reference  proper  motions 
(not  positions)  when  these  are  available.   Then  we  have,  in 
the  same  notation  as  above,  a  set  of  observation  equations 
for  the  m-th  plate, 


^m^^'y  Pm^  =  ^^o-^^'^m^  "^   ^  <3) 


where  now  \i       is  known  for  some  of  the  stars.   Again  a  least- 
squares  adjustment  for  the  p's  is  made  for  each  of  the  M 
plates.   Then,  for  each  of  the  remaining  stars  on  the  m-th 
plate. 


^m  =  ^^m^^'^'P^-^o^/'^m'  (^> 


m 
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In  order  to  take  full  advantage  of  the  material,  one  should 
then  treat  these  derived  y's  as  known  values  in  (3)  and 
iterate . 

This  procedure  is  still  subject  to  the  restriction 
that  all  plates  be  taken  with  the  same  telescope  and  with 
the  same  plate  centers,  and  that  there  be  a  sufficient  num- 
ber of  stars  in  the  field  with  known  proper  motions  to 
determine  the  six  adjustment  constants. 

Classical  Plate  Reduction 

We  can  depart  from  the  requirement  that  all  plates  be 
taken  with  the  same  telescope  with  stable  optical  proper- 
ties and  with  the  same  plate  centers  (termed  together  "pro- 
jection conditions")  by  referring  the  measured  coordinates 
to  a  system  of  reference  star  positions  instead  of  to  one 
another.   Let  the  spherical  coordinates  of  the  stars  pro- 
jected onto  a  plane  approximately  tangent  to  the  plate  cen- 
ter, the  conventional  "standard  coordinates,"  be  (C»n). 
Then  the  observation  equation  relating  these  to  the  mea- 
sured coordinates  is,  following  the  notation  used  above. 


f    (?+M    t    ,    n+y    t    {p}    )    =   X  +Ax 
m  xm  ymm  o 


(5) 


where  f   is  in  general  some  polynomial  in  E, ,     n,  and 
magnitude.   The  f   is  now  complicated  by  the  appearance  of 
terms  expressing  tilt,  coma,  radial  distortion,  etc.   (See, 
for  example,  Konig,  1962  and  Eichhorn,  1963  for  general 
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discussions.)   (Here,  we  have  assumed  that  C  and  n  vary 
linearly  with  time;  this  is  an  approximation  valid  for  the 
small  changes  involved.) 

These  equations  are  solved  for  the  {p}   by  separate 
least-squares  adjustments  for  each  plate.   Then  for  each  ol" 
the  non-reference  stars  we  have,  for  each  plate  the  vector 
equation  (where  E,    now  stands  for  the  vector  (C,ri)), 


m     m 


from  which  the  proper  motions  are  computed  by  solving  for 
each  star  the  set  of  ra  equations 


C   +  yt   =  C 
o     m    m 


There  are  two  advantages  to  the  plate  reduction  method 
over  the  differential  methods.   First,  it  does  not  require 
that  the  plates  have  been  taken  with  the  same  telescope  and 
plate  centers.   Thls_  in  itself  is  the  decisive  factor  when 
the  plate  material  is  unavoidably  heterogeneous,  as  in  the 
present  program.   Second,  tbs  plate  reduction  yields  actual 
star  positions,  in  addition  to  the  proper  motions. 

The  disadvantage,  as  stated  above,  lies  in  the  neces- 
sity of  including  higher  terms  in  the  observation  equation. 
The  presence  of  these  terms  necessarily  increases  the  total 
error  due  to  all  the  parameters.   Consequently,  star  posi- 
tions which  are  computed  from  them  contain  large  systematic 
errors,  especially  toward  the  edge  of  the  plate  (Eichhorn 
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&  Williams,  1963)  and  also  in  the  spaces  between  widely 
separated  reference  stars.   This  type  of  error  ("parameter 
error")  is  not  revealed  by  the  adjustment  residuals  of  the 
reference  stars,  but  only  by  comparing  the  positions  of 
field  stars  obtained  from  several  plates  with  an  indepen- 
dent position,  or  with  each  other.   The  more  closely  dupli- 
cated are  the  circumstances  of  the  plates,  the  less  apparent 
are  these  parameter  errors. 

Another  weakness  of  the  conventional  plate  reduction 
procedure,  implicit  in  equation  (5),  is  the  assumption  that 
only  the  plate  measurements  are  subject  to  error.   In  fact, 
however,  the  reference  star  positions  are  subject  to  error, 
too.   For  example,  a  position  error  of  '.'3  (not  at  all 
atypical  for  the  Yale  Catalogue)  corresponds  at  the  focal 
plane  of  a  typical  "parallax"  telescope,  to  15ii--which  is 
far  larger  than  the  expected  measuring  error. 

One  might  therefore  be  tempted  to  treat  each  plate  as 
a  separate  adjustment.   Approximately,  this  can  be  accom- 
plished by  assigning  the  x  and  C  relative  weights  according 

to  some  assumed  a  ,  a^.   All  that  this  accomplishes,  how- 

x'   C 

ever,  is  to  divide  each  residual  into  a  "star"  component  and 
a  "measurement"  component  in  the  ratio  o    /a    ,  effecting  a 
mainly  spurious  reduction  in  the  measurement  residuals. 
Instead,  it  is  much  more  important  to  enforce  the  condition 
that  the  star  position  correction  for  a  given  star  be  a 
single  value,  not  a  set  of  values  for  each  plate  on  which 
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that  star  appears.   It  is  the  rigorous  expression  of  this 
condition  which  prompted  the  development  of  the  plate  over- 
lap method  treated  below. 

Among  the  proper  motion  investigations  employing  the 
classical  plate  reduction  method  there  are  two  of  the  Orion 

Nebula  cluster:   Meurers  and  Sandmann  (1963)  and  (for  some 

of  the  material)  Parenago  (1954). 

Differential  Overlap 

We  pause  to  remark  here  that  the  overlap  condition  can 
be  applied  in  a  differential  formulation  as  well. 

In  the  notation  used  above,  we  have  for  the  1st  star 
on  the  first  and  m-th  plates 


^l^^l'^l'^-P^l^ 


t,Ay^ 


^m^^r^l'^Pm^)  +  "^m^ 


t  Au, 
m  1 


\  +   t^q^  +   Ax     (6a) 


=  X  +  t  d,  +   Ax 
X    ml 


0    +  t^Ap^    (6b) 


and    for    the    i-th    star    on    the    m-th    plati 


^m^^'^i'^Py 


mi  1  m   1  1 


t  Ay.    = 
m     1 


+  t  Ay. 
m     1 


Again  we  compute  the  {p}^  and  Ap^  for  the  reference  stars  by 
a  least-squares  adjustment  and  by  successive  iterations,  the 
y's  of  the  field  stars.   In  forming  the  normal  equations, 
the  conditions  (6b)  are  weighted  with  respect  to  the 
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2      2 
observations  (6a)  by  a  factor  a  (x)/a  (y)  chosen  a  priori. 

The  matrix  of  the  normal  equations  can  be  reduced  to 
dimension  equal  to  the  number  of  Ap ,  as  described  in  the 
section  "plate  overlap."   The  Ap  will  be  distributed  nor- 
mally with  mean  zero  and  dispersion  =  a(y)/a(x)  times  the 
dispersion  in  Ax. 

Of  all  the  formulations,  this  should  give  the  most 
accurate  proper  motions;  but  to  be  applicable,  the  plate 
material  must  be  homogeneous  and  there  must  also  be  an  ade- 
quate number  of  stars  with  known  proper  motion  within  the 
field  of  view.   These  conditions  are  not  satisfied  in  the 
present  program. 

To  the  author's  knowledge  this  formulation  has  never 
before  appeared  in  the  literature. 

Plate  Overlap  Method 

The  Observation  Equations 

We  now  write  the  observation  equations  in  their  most 
general  form.   In  order  to  illustrate  the  pattern  more 
clearly,  we  give  here  the  equations  for  a  reference  star 
(the  first  star)  and  a  field  star  (the  i-th  star)  as  they 
occur  on  the  first  and  ra-th  plates. 


•J--  > 
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f^d^+t^fi^Av)^  +  ^i^^i+^itiAy^  =(x^+Ax^)F^    (7a) 

f   (|^+t^11^,{p}^  +  F^AC^+t^AMj^=(x>Ax^)F^    (7b) 

f   (L+t   0,    ,{p}     +  F  A5,+F   t  Ap  =(x,+Ax.)F^ 

mlmlxm         mlmml  11m 

f   il.+t  a.    ,{p}  F  AC.+t  Am.=(x.+Ax,)F 

ml     m  ix'   "^m  mimiilm 

+       ACj^  =  O+A?^  (7c) 

"tAy  =  0+tAy^  (7d) 

where  for  brevity  f(?,{p})  stands  for  two  functions,  i.e., 
is  a  1  X  2  matrix,  and  x,)j,?,{p}  stand  for  the  vectors 

(x.y).  ^^x'^y^'  (?.ri)>  (P]^.P2---)  respectively.   F^  is  the 
reciprocal  of  the  assumed  focal  length  in  the  same  units  as 
the  X,  of  the  telescope  with  which  the  m-th  plate  was 
taken.   (While  the  focal  lengths  are,  in  effect,  among  the 
parameters  being  adjusted,  only  very  approximate  values  are 
needed  for  this  role  and  they  need  not  be  revised  in  sub- 
sequent iterations.)   The  E,    and  p  are  in  units  of  radians 
and  radians/yr.,  respectively.   Carets  denote  the  values 
used  in  the  initial  approximation. 

The  tm  are  the  plate  epochs  in  years  now  measured  from 
the  reference  catalogue  epoch,  in  the  case  of  reference 
stars,  and  from  some  arbitrary  epoch  in  the  case  of  field 
stars.   We  have  incorporated  the  approximation  x(t)  = 

-r~  Mt  =  Fy  t,  fy  t  into  the  terms  involving  y. 
9e         X      y  ^ 

In  the  right  hand  colum  appear  the  observed  quanti- 
ties.  In  the  equations  (7c,  7d)  representing  the  reference 


'iW«f-'.- 
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positions  and  motions,  the  observed  quantities  are  the  de- 
partures of  the  initial  values  of  the  reference  position 
or  motions  from  the  assumed  values;  these  are  equal  to 
zero.   The  A  appearing  in  the  right  column  symbolize  the 
errors  of  the  observations;  on  the  left  side  the  A  are 
unknowns.   The  "Aia"  of  the  field  stars  are,  ordinarily, 
the  y  themselves.   Of  course  there  must  be  at  least  two  ob- 
servations at  sufficiently  separated  t   for  each  field  star 

m 

whose  y  is  to  be  determined.   If  this  condition  is  not  met 
for  a  particular  star,  we  merely  eliminate  its  u    from  the 
equations . 

In  order  to  handle  the  entire  set  of  equations  by  the 
conventional  least-squares  algorithm,  the  equations  are 
scaled  so  that  the  scaled  A  are  members  of  one  normal  dis- 
tribution with  a  single  standard  deviation  o  .   Here  this 

"  u 

standard  deviation  is  measured  in  radians,  hence  the  factor 
F  in  the  equations  (7a).   The  proper  motion  equations  (7c) 
are  each  scaled  by  t,  a  time  value  such  that  to   =  a  . 
Finally,  the  non-uniformity  of  the  a    of  the  various  cata- 
logues and  of  the  plate  measurements  is  accounted  for  in 
the  usual  manner  by  multiplying  each  equation  by  the  appro- 
priate value  0  /a(  =  /w). 
u 

The  Normal  Equations 

The  matrix  of  the  resulting  system  of  normal  equations 
is  a  single  large  array  of  dimension  equal  to  the  total 
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number  of  unknowns--rather  than,  as  in  the  classical  re- 
duction, a  set  of  matrices  for  each  plate  of  dimension  equal 
only  to  the  number  of  plate  constants  per  plate.   This  total 
number  of  unknowns  is  very  large  for  the  number  of  plates 
and  stars  involved  in  a  typical  proper  motion  study;  the 
program  which  is  the  subject  of  this  paper,  involving  88 
plates  and  130  stars,  leads  to  some  2,000  unknowns.   Clearly 
this  is  an  unmanageably  large  dimension.   We  therefore  have 
taken  steps  to  reduce  the  size  of  the  largest  matrix  which 
must  be  inverted. 

First,  we  may  always  align  the  plates  in  the  measuring 
machine  so  that  the  instrumental  x-axis  is  nearly  parallel 
to  the  5  axis.   Then 

— =  TVrAf  -4-  ^■^,      \-         ^ 


Thus  the  matrix  E  consists  of  elements  of  dimension  2x2 
rather  than  4x4. 

Second,  we  assume  the  plate  constants  p  appearing  in 
the  equation  f(5,n>{p})  =  x  to  be  independent  of  those  in 
f'(5»n,{p'})  -  y.   This  means  that  the  normal  equations 
involving  the  unknowns  p.  A?,  Ay   are  decoupled  from  those 
in  p'.  An,  An   leaving  two  independent  sets  of  mormals  each 
of  half  the  rank  of  the  original. 

The  system  of  normal  equations  may  be  represented 
schematically  as  follows.   Let  p   represent  submatrices  due 
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only  to  coefficients  of  p  (plate  constants  of  m-th  plate)  ; 
E  represent  the  matrix  due  only  to  coefficients  of  star  un- 
knowns (AC, Ay):   this  is  a  block  diagonal  of  2  x  2  blocks; 
C,  the  matrix  of  cross-terms,  and  V^(i)  the  vector  of  plate 
constant  coefficients  for  the  i-th  star,  on  the  m-th  plate. 
On  the  right  side,  the  vector  L  consists  of  the  homologous 
vectors  of  the  classical  reduction  normals,  strung  together, 
while  the  elements  of  x  are  given  by  x^  =  T  W^  F^  X^^ . 


?1 

n 

?2 

c 

p 

0         p 

m 

C^ 

n 
u 

0  ^-q 

In  an  earlier  application  of  the  overlap  method,  Googe 

(1967)  showed  that  the  system  of  normals  can  be  reduced  in 

dimension  to  only  the  total  number  of  plate  constants  by 

"folding,"  that  is,  by  eliminating  E  and  filling  up  the 

entire  block  containing  the  P  .   (In  that  application  E  was 

m 

a  long  diagonal  matrix  of  single  elements.)   For  the  task 
at  hand,  we  instead  eliminate  the  P,  folding  the  matrix 
into  E  and  producing  a  solid  block  of  dimension  equal  only 
to  n  X  n,  n  being  the  number  of  stars.   Since  n  =  130  (as 
opposed  to  about  ten  times  that  number  of  plate  constants) , 
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we  are  left  with  a  matrix  that  can  be  easily  handled  and 
Inverted  with  present  computer  capacities. 

These  reduced  matrices,  S  and  T,  can  be  formed  in 
place  without  the  need  of  auxiliary  matrices.   They  are 
simply 


S  =  E3_  -  C^P-^C^^ 


T  =  E^  -  C^P-'C^^ 


where  the  subscripts  merely  refer  to  the  homologous  counter- 
parts in  the  matrices  of  the  original  normal  equations. 

It  can  be  seen  that  the  matrix  C  is  composed  of  strings 
of  vectors  such  that  the  i-th  row  of  C  is 

V,   (i),  V^   (i) ,  ...  V    (i) .   Thus,  the  elements  of  S 
l,p       2 ,p  m,p 

(and  similarly  of  T)  can  be  written 


M 


S.  .  =  w.(6.  .)  -  ^  V  (i)P"^V^(j) 


m 


m 


m  m 


with  a  similar  equation  for  T,  ..   The  Kroenecker  delta 

denotes  the  original  diagonal,  V  (i),  is,  again,  the  vector 

of  plate  constant  coefficients  for  the  m-th  plate  and  1-th 

star,  and  w.  are  the  star  weights.   For  simplicity  we  have 

ignored  the  proper  motion  components  of  S.   Thus  the  entire 

array  can  be  formed  from  one  plate  at  a  time  and  built  up 

In  layers.   No  more  data  need  be  stored  at  a  time  than  the 

measurements  for  a  single  plate,  the  inverse  P    of  the 

m 
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matrix  of  normal  equations  generated  by  the  single  plate  m, 
and  the  arrays  S  and  T  with  the  updated  sums  as  their  ele- 
ments.  Only  the  upper  halves  of  S  and  T  need  be  maintained, 
since  they  are  symmetric.   Another  convenient  feature  of 
the  formulation  is  that  the  data  generated  by  new  plates 
may  be  added  as  they  become  available  without  the  necessity 
of  completely  recomputing  the  already  existing  S  and  T. 


SECTION  V 
THE  COMPUTATIONS 


Partition  of  the  Plate  Material  and  Treatment 
of  the  Proper  Motion 

The  plate  material  naturally  divides  itself  into  three 
groups:   an  early  (^1900  group  (A)  of  wide-field  photo- 
graphy; a  later  (=1905-1927)  group  (B)  of  long-focus  photo- 
graphy; and  the  recent  group  (1968-1974)  of  long-focus 
photography.   We  have  taken  advantage  of  this  natural  group- 
ing to  effect  a  certain  measure  of  simplification.   The 
procedure  adopted  is  (1)  Reduction  of  the  recent  group  in 
a  single  overlap,  without  regard  to  the  proper  motion  un- 
knowns, using  as  initial  values  for  the  field  stars  the 
positions  calculated  from  plate  constants  obtained  for  a 
subset  of  these  plates  reduced  using  reference  stars  only; 
(2)  Reduction  of  the  early  group  as  a  single  overlap,  also 
without  regard  to  the  unknown  p;  (3)  Calculation  of  a 
"secondary  catalogue"  of  a,  6,  y^  ,  \i ^    and  means  epochs  for 
the  field  stars,  using  the  results  from  steps  (1)  and  (2); 
(4)  Reduction  of  group  (B)  using  this  "catalogue"  of  field 
stars  as  reference  material  of  lower  weight,  in  addition  to 
the  original  catalogue  of  positions  and  motions;  (5)  Cal- 
culation of  new  proper  motions  of  the  field  stars  from  the 
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solutions  of  steps  (1)  and  (4) ;  and  (6)  interaction  of 
steps  (4)  and  (5),  resulting  in  the  final  values  of  the 
proper  motions.   This  procedure  is  justified  by  the  narrow- 
ness of  the  range  of  epochs  in  groups  (A)  and  (C),  while 
the  22-year  spread  of  group  (B)  is  adequately  compensated 
for  by  the  use  of  the  approximate  proper  motions  for  its 
field  stars . 

The  role  played  by  the  Astrographic  Catalogue  and 
Z6-Se  plates  is  in  effect  to  provide  starting  values,  well 
determined  with  regard  to  scale  and  introduced  systematic 
errors,  for  the  positions  and  motions  of  the  field  stars 
appearing  on  the  Yerkes  and  McCorraick  p lates-- these  latter 
plates  containing  few  reference  stars.   The  first-named 
plates  are  well  suited  for  this  role  because  they  cover  a 
wide  area  including  many  reference  stars.   (See  Table  2.) 
The  relatively  small  scale  of  these  plates  is  no  great  dis- 
advantage here  since  high  internal  accuracy  is  not  essen- 
tial for  the  starting  values.   (We  recall  that  these  ini- 
tial values  are  treated  as  reference  material,  not  merely 
as  starting  values.) 

In  contrast,  the  late-epoch  photography  (group  (C)), 
being  overlapped,  contains  a  sufficient  number  of  reference 
stars  to  allow  a  straightforward  reduction  without  the  need 
for  a  "secondary"  catalogue  (but  see  below,  "Stability  of 
the  Solution") . 
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Behavior  of  the  Solution 
We  need  not  dwell  on  the  routine  mechanics  of  the 
calculations.   We  are  concerned  here  with  properties  of  the 
solution  which  are  only  revealed  when  the  calculations  are 
actually  made.   To  the  extent  that  such  properties  are  not 
completely  predictable  in  advance,  we  may  say  that,  in 
effect,  they  reflect  the  stability  of  the  solution.   It 
ought  to  be  remarked  at  this  point  that  there  has  not  pre- 
viously appeared  a  plate  overlap  reduction  of  the  general- 
ity which  obtains  here,  wherein  the  number  of  unknowns 
greatly  exceeds  the  number  of  reference  stars  as  is  typ- 
ically the  case  in  proper  motion  studies.   The  following 
aspects  of  the  behavior  require  special  mention. 


Iterations 

A  solution  of  the  normal  equation  yields  a  set  of  star 
position  corrections,  and  plate  constants  obtained  from  the 
initial  positions.   One  therefore  replaces  the  initial  star 
positions  for  the  revised  values  and  iterates  the  solution 
until  all  the  star  position  corrections  are  less  than  some 
value.   It  was  found  that  convergence  is  quite  slow,  re- 
quiring about  seven  iterations  for  convergence  to  within 
V02  for  all  the  star  positions.   Convergence  is  more  rapid 
when  fewer  field  stars  are  included  as  unknowns.   If  no 
field  stars  are  included,  the  solution  convergences  with 
only  two  iterations. 
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Influence  of  Weighting 

It  was  found  that  the  solution  is  rather  sensitive  to 
relative  weights  assigned  to  the  several  source  catalogues. 
If  too  high  a  weight  is  assigned  to  a  star,  the  residuals 
for  that  star  from  the  various  plates  on  which  it  appears 
reveal  a  consistent  error.   Thus  it  was  found  that  the 
accuracy  of  positions  and  motions  from  the  FK4  and  N30  is 
lower  than  the  quoted  mean  errors  for  those  catalogues. 

If  too  low  a  weight  is  assigned  to  a  star,  there  is 
little  effect  if  the  catalogue  value  is  very  close  to  the 
true  value.   If  the  initial  position  is  incorrect  by,  say, 
0V4,  however,  the  reduction  will  over-correct  its  position. 
Moreover,  it  can  be  argued  intuitively  that  if  a  reference 
star  occurs  in  a  relatively  isolated  region,  the  consequence 
of  assigning  its  position  too  low  a  weight  is  to  allow  its 
position  to  adjust  by  too  far,  without  the  over-adjustment 
being  evident  in  the  residuals.   This  is  expected  because 
the  plate  constants  are  freer  to  absorb  an  error  in  a  star 
position  if  there  are  no  surrounding  stars  to  help  con- 
strain the  values  of  the  plate  constants  evaluated  in  that 
region.   Indeed,  it  was  found  by  experiment  that  when  the 
position  of  one  reference  star,  15'  from  its  nearest  neigh- 
bor, was  allowed  to  vary  over  a  range  of  O'.'S,  the  residuals 
for  that  star  varied  by  no  more  than  Vll   The  implication 
to  be  drawn  from  this  experience  is  that  it  is  preferable 
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to  overestimate  rather  than  underestimate  the  accuracy  of 
the  catalogue  positions. 

The  weight  assigned  to  the  plate  measurements  rela- 
tive to  the  catalogue  positions  also  were  found  to  exert  a 
sensitive  influence  on  the  solution.   We  do  not  know  in  ad- 
vance exactly  the  mean  error  to  be  assigned  to  the  plate 
measurements,  especially  those  taken  from  the  literature. 
Within  reasonable  limits  (4vi  to  8p),  our  choice  of  the  mea- 
surement mean  errors  is  guided  by  the  conjecture:   given 
several  sets  of  observations  each  characterized  by  an  inde- 
pendent error  distribution  with  some  0,  that  set  of  a  which 
minimizes  the  variance  of  the  adjustment  made  on  the  obser- 
vations is  the  best  obtainable  estimate  of  the  a. 

Weighting  of  Field  Stars 

Since,  by  definition,  the  field  stars  have  no  a  priori 
positions,  the  initial  values  assigned  to  them  have  zero 
weight;  that  is,  the  observation  equations  (7c) 


A5^/w^  =  AC./w^ 


corresponding  to  them  reduce  to 


0  =  0 


It  was  found  that,  in  practice,  some  non-zero  weight 
must  in  fact  be  assigned  to  the  field  star  positions;  other- 
wise the  solution  is  so  ill-conditioned  that  successive 
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iterations  rapidly  diverge.   This  ill-conditioning  worsens 
as  the  number  of  field  stars  surpasses  the  number  of  refer- 
ence stars,  but  improves  with  improvement  in  the  initial 
values . 

We  have  therefore  been  compelled  to  assign  a  weight  to 
the  initial  values  of  the  field  star  positions  correspond- 
ing to  1/10  that  of  the  Yale  positions  at  the  same  epoch. 

Residuals 

There  is  a  residual  for  each  star  on  each  plate.   The 
residuals  obtained  from  the  final  iterations  have  the  fol- 
lowing rms  errors,  calculated  plate  by  plate,  shown  here 
for  each  telescope.   The  values  obtained  from  all  the 
plates  except  a  few  anomalous  cases  are  included  in  the 
ranges  given. 


TABLE  3 
RMS  ERRORS 


Telescope  rms 

Yerkes  V09-V06 

Allegheny  .09-. 07 

McCormick  .08-. 06 

USF  .09-. 04 

San  Fernando  .40-. 25 

Zo-Se  .35-. 15 
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We  also  compute  for  each  star  the  average  of  its  resid- 
uals  for  each  plate.   This  mean  residual,  added  to  the 
calculated  position  correction  and  the  initial  position, 
represents  the  position  for  that  star  as  calculated  from 
the  final  values  of  all  the  plate  constants.   It  is  this 
value  which  is  taken  as  the  final  calculated  position  of 
the  star  at  the  corresponding  mean  epoch  of  the  plates  on 
which  it  appears. 


The  Computer  Program 
The  computer  program  with  which  these  calculations 
were  carried  out  is  written  in  Fortran  IV.   It  is  con- 
tained in  only  ahout  800  statements  and  is  conveniently 
modular  in  structure.   It  utilizes  two  tracks  of  working 
disk,  but  requires  no  tape  storage.   Taking  advantage  of 
the  cumulative  character  of  the  matrices  of  normal  equa- 
tions, the  program  accepts  an  indefinite  number  of  plates, 
with  an  arbitrary  (up  to  the  maximum  value  of  n)  number 
of  stars  per  plate,  in  arbirary  order.   The  number  n  of 
stars  itself  is  completely  arbitrary,  within  the  limita- 
tions of  storage  capacity.   With  core  storages  of  180K, 
240K,  and  300K,  the  corresponding  maximum  values  of  n  are 
80,  112,  and  135.   Were  the  proper  motion  explicitly  in- 
cluded, the  corresponding  values  of  n  would  be  1/2  of  these 
values.   The  catalogue  of  reference  stars  is  treated  as  an 
interchangeable  data  set.   Copies  of  the  program  on  cards 
are  available  to  interested  persons  on  request. 
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The  Plate  Models 

By  plate  model  we  mean  the  polynomial  terms  required 
to  transform  the  Standard  Coordinates  of  the  stars  into 
their  (x,y)  coordinates  on  the  plate.   Each  telescope  can 
be  expected,  in  general,  to  require  its  own  model;  finding 
the  proper  model  for  a  given  telescope  is  a  matter  of 
experiment . 

For  all  plates,  the  terms  whose  coefficients  are  1, 

2  2 

C»  n  >  C  >  Cl  (i"  the  x-equation)  and  1,  E,  ,    n  .  n  >  ?n  (in  the 

y-equation)  are  included,  since  these  terms  always  arise 

from  the  projection  geometry.   The  plate  models  finally 

adopted  for  each  telescope  include,  in  addition  to  the 

above  terms,  the  following. 


TABLE  4 
PLATE  MODELS 


Terms  (and  origin) 

2     2       2   2 
m         xm ,  ym      xr  ,  yr       y  ,x 

Telescope       (guiding)     (coma)     (radidist)      (?) 

USF  X  X 

Allegheny  X  X 

Yerkes  XX  X 

McCormick  X  X 

Z5-Se  X  X 

San  Fernando  X  X 


i't^^''-/ 
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No  determination  of  terras  involving  color  could  be 
made,  even  though  such  terms  are  undoubtedly  present, 
because  the  distribution  of  colors  is  Loo  strongly  cor- 
related with  magnitude  and  position.   Most  of  the  program 
stars,  while  not  members  of  the  Trapezium  cluster  proper, 
are  members  of  the  Orion  aggregate;  consequently,  the  bluer 
stars  are  also  brighter.   Moreover,  nearly  all  the  red 
stars  are  concentrated  in  the  near  vicinity  of  the 
Trapezium.   An  attempt  to  model  color  terms  in  these  cir- 
cumstances could  result  in  grave,  undetectable  systematic 
errors . 

It  was  found  that  the  USF  telescope  shows  neither  sig- 
nificant coma  nor  radial  distortion;  some  plates  show  a 
linear  magnitude  term  in  x  or  y.   The  Allegheny  telescope 
(Thaw  refractor)  also  shows  negligible  coma,  but  does  ex- 
hibit some  radial  distortion;  the  earlier  plates  tend  to 
exhibit  larger  magnitude  terms  than  do  the  recent  plates. 

The  Zo-Se  plates  show  surprisingly  small  higher-order 
terms . 

The  McCormick  plates  show  both  linear  magnitude  terms 
and  large  and  variable  coma,  apparently  amounting  to  about 
Ip/cm-mag-  But  the  determination  of  the  coma  terms  for 
these  plates  is  complicated  by  several  factors:  the  stars 
appearing  on  them,  especially  away  from  the  center,  have  a 
narrow  range  of  magnitude,  and  they  were  taken  without  an 
objective  grating.   Moreover  there  is  a  strong  correlation 
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between  the  coma  term  and  the  magnitude  term  in  declina- 
tion.  Rather  than  risk  the  introduction  of  unnecessary 
parameter  error,  we  finally  pre-correct ed  the  measurements 
for  the  effect  of  coma,  using  for  the  purpose  the  coma 
value  found  earlier  by  Eichhorn  et  al   (1970)  from  photo- 
graphy which  was  better  suited  than  ours  for  its 
de terrain at  ion . 


SECTION  VI 
RESULTS 

The  derived  positions  and  proper  motions,  given  for 
orientation  1950,  are  listed  in  Table  5.   The  right  ascen- 
sions and  declinations  are  given  in  decimals  of  degrees  as 
the  best  compromise  between  convenience  of  computation  and 
of  identification.   All  the  positions  are  given  for  epoch 
1974.0,  this  being  close  to  the  actual  epoch  of  the  recent 
plates.   The  weighted  mean  epoch  of  the  recent  plates  is  also 
given  for  each  star.   The  proper  motion  accuracies  are  cal- 
culated from 

a2p  =  ((j2    +    a2)/At   , 
1     2 

where  o^  and  a^    are  the  mean  errors  of  position  at  the  mean 
early  and  mean  recent  epochs.   These  in  turn  were  calculated 
simply  from  the  average  mean  errors  a    of  the  n  plates  on 
which  the  star  appeared  at  each  epoch  according  to 

Absolute  Proper  Motion  of  the  Trapezium 
The  average  proper  motion  of  the  Trapezium  Cluster  is 
li^  =  VOlO/yr,  M^    =  VOOl/yr 
Part  of  this  motion  is  due  to  secular  parallax,  which  amounts 
to  (assuming  a  distance  of  500  parsecs) 

y„  =  +V001/yr,  Mg  =  -''OOS/yr 
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The  galactic  rotation  terra  in  proper  motion  in  t'ae 
direction  of  the  Orion  Nebula  amounts  to  less  than  V0005/yr 

Velocity  Dispersion 

The  most  remarkable  result  emerging  from  the  proper 
motions  is  the  very  small  value  of  the  velocity  dispersion. 
This  is  given  in  Table  6  for  various  sub-groups  of  stars. 
Here  "a"  denotes  dispersion,  not  mean  error. 

Table  6 


No.  of 
Group  Radius     Members  a       C'/yr)    a       C'/yr) 


Trapezium 

el  and  9^ 

Trapezium  Cluster 
Nebula  Cluster 


(a):   excludes  stars  assumed  to  be  foreground  objects. 

The  dispersion  in  proper  motions  for  the  inner  members 
of  the  cluster  is  so  small  as  to  call  into  question  whether 
any  true  motion  is  detectable  at  all.   In  any  case  it  sets 
a  realistic  upper  limit  to  the  true  proper  motion  errors, 
since  random  errors  will  act  to  increase  the  observed  dis- 
persion.  It  is  believed  that  these  are  the  smallest  disper- 
sions ever  obtained  for  a  group  of  stars. 

Translating  the  proper  motions  into  tangential  velocities 
at  the  assumed  distance  of  500  parsecs,  we  obtain  a  velocity 


% 

1^6 

10" 

4 

.0010 

.0009 

2' 

9 

.0013 

.0012 

6' 

43 

.0015 

.0014 

20' 

93^ 

a) 

.0036 

.0026 
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dispersion  of  3h    km/sec  (in  one  coordinate)  for  the 
Trapezium  Cluster,  Ih   km/sec  for  the  larger  Nebula  Clus- 
ter.  (It  is  worth  recalling  that  the  velocity  dispersion 
of  the  Pleiades,  which  is  less  than  1/3  the  distance  of 
Orion,  is  about  l^j  km/sec.)  We  suspect  that  part  of  the 
difference  in  the  two  figures  is  spurious,  due  to  the  increase 
of  error  as  one  moves  away  from  the  region  near  the  plate 
centers  and  greatest  concentration  of  stars.   It  is  interes- 
ting that  Strand  obtained  an  equally  small  proper  motion  dis- 
persion for  the  larger  Nebula  Cluster,  but  a  much  larger 
dispersion  for  the  Trapezium  and  the  Trapezium  Cluster 
('.'008/yr.,  and  V004/yr.,  respectively).  This  is  probably  due 
to  the  fact  that  on  the  Yerkes  plates  the  nebula  and  the 
brighter  stars  are  overexposed. 

Contraction 
For  each  proper  motion  let  the  radial  and  tangential  com- 
ponents be  ]i      and  y  ,  the  radial  distance  R  being  reckoned 
R       T 

from  the  Trapezium.   Let  the  coordinates  of  the  star  measured 
parallel  to  a  and  6,  with  the  Trapezium  at  the  origin  be  A,  D. 
A  linear  expansion  (or  contraction)  can  be  represented  by  a 
coefficient  k  in  units  ("/year) /R(degrees) .   The  expansion's 
contribution  to  \i-d    is  an  amount  u',  such  that 

R     R     R 
Then  in  order  to  find  k,  we  solve  by  least  squares  the 

equations 

R^k  +  (A/R)^y^  +  (D/R)^7j,  =  (y^)i  +  (AyR)i 
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where  'ii^,    vTa  Is  the  mean  proper  motion  of  the  group.   If 
y  ,  y   are  not  known,  a  priori  and  if  the  distribution  of 
stars  is  non-uniform,  they  also  must  be  treated  as  unknowns. 

The  value  of  k  determined  from  our  proper  motion  in  the 
Nebula  Cluster,  excluding  the  high  proper  motion  star  is 

k  =  (V009  +  .002/yr) /degree   ; 
that  is,  a  contraction  of  the  cluster  amounting  to  V003/year 
at  its  edge.   This  result  was  quite  unexpected  and  prompted 
a  thorough  search  for  sources  of  error  that  might  have  pro- 
duced it  as  an  artifact.   As  a  weak  check  on  the  result  we  can 
solve  for  contraction  in  a  and  in  5  separately  (such  "con- 
tractions" having  no  physical  significance  considered 
separately);  it  is  found  that  the  two  coefficients  are  equal 
within  one  standard  deviation.   It  must  be  borne  in  mind  that 
the  uncertainty  of  this  result  is  greater  than  indicated  by 
its  formal  error.   This  is  because  the  contraction  coefficient 
depends  critically  upon  the  proper  motions  of  the  stars  at 
the  edge  of  the  cluster  (which  in  itself  is  indeed  reflected 
in  the  formal  uncertainty  of  the  result),  while  it  is  the  posi- 
tions and  hence  proper  motions  of  precisely  these  stars  which 
are  most  afflicted  (as  pointed  out  above,  p.  40)  and  to  an 
uncertain  extent  by  the  plate  constants  errors.   In  any  event, 
let  us  accept  our  result  and  examine  its  consequences. 

Translated  into  linear  velocity  at  the  edge  of  the 
Gerola-Sofia  cloud,  the  R=4  parsecs,  the  contraction  amounts 
to  10  km/sec.   This  is  in  surprising  agreement  with  the  value 
12  km/sec  obtained  by  Gerola  and  Sofia  for  their  model  of  the 


*!ikf^' 
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cloud  and  lends  strong  support  to  their  interpretation  of  the 
cloud  molecular  line  widths. 

Quite  independent  of  the  contraction  determination  is 
the  evidence  of  the  velocity  dispersion  in  the  Trapezium 
Cluster.   Since  that  result  argues  against  the  existence  of 
turbulent  eddies  with  velocities  much  higher  than  3  km/sec, 
which  seem  to  be  required  if  the  cloud  is  not  contracting, 
we  have  a  double-barrelled  support  of  the  view  that  the  Orion 
Nebula  is  relatively  quiescent,  contracting  under  gravity, 
and  has  not  yet  entered  its  most  productive  era  of  star 
formation . 

If  the  cluster  is  not  expanding,  the  "expansion  ages" 
assigned  to  it  in  the  post  are  invalid,  and  stars  already 
formed  in  the  cluster  cannot  be  assigned  a  single  age  on  the 
basis  of  kinematic  considerations.   The  contraction  removes 
any  discrepancy  between  the  kinematic  and  other  age  determi- 
nations, and  allows  for  a  continuous  range  of  ages  in  the 
member  stars. 


68 


oor--iriunvo<-in'-H<f<f\£)u^inr-^in«a-Ln-;fLr»\CJ<rinvD<-<r<r<!-<riri 

•-HOCCOOO'-HOCOOOOOCOOOOOOOOOOOOC 

c 
o 


u 
o 

C 

ta 

r. 


r->oooi— looc^lo^c^lvoc^Joos3■-d• 
•-^r~m^ooc^Jl^)c^^u^l— iiHinm 


v3-o-sj-ccLn<rcsicNcM 
inmincvjinmoocviiH 


cTi  c»  o)  o-v  m 
«*  i-i  r-N  ^-  cvj 


o^  cjN  00 
vn  s"  vo 


crv  CM  o  cri  00 
vo  r^  I —  vo  vo 


O^CN««3-CSOC^iO 


c^i  O  ro  n  cy 
1^  rv  r>  r^  vc 


en  c  cr>  o^  fO 
r-~  r^  kC  \o  r~ 


O 


U 


ocvioovo-3-vor^o><T(vovor~-fOvo<riri<ras 
cyiiTi'-Hc^omoOin'-tcviiHiHoiooooco 

vrOOOOOOO—iOOOOOOrHOO 

ocoooocooooooooooc: 


iHrHsrOOCMtHiHinOOrH 

coooooooooo 
ocoooooooco 


I  I 


I 


! 


I  I 


I 


I  I 


I 


I  I  I 


I  I 


cvj<-cy^-^mi~'<-<Lnoomvoc>ivovoc^-j-ocovocnrHcovovcr^vor~o>iH 


r-HO'-<r-CN)OCO'-lvDfOCTiCV)LriCN\OiHvOrH 
OCMOOOr-io<-HOOOFHOOOCOO 

oooooooooooooooooo 


vOl-IC^00a^^~mo^C•.  1— I 
OOOOOOOOOiH 

ooooooocoo 


t    I 


I 


I    I    I    I 


I    I    I 


I    I 


I     r    I 


I    I    I 


U1 


H 

o 

s 

w 
p^ 
o 

W  Pri 

h-l  Pi 

<!  O 

E-l    |S5 


CO 

5?; 
o 

M 

H 
M 

o 

P:< 


incoor~-c-iooo(Jcr>mnvocMOioirirMO\cncNivooooo>inms;j-o\<i- 

i-iN    lO    rri    »-»•    ir»    rf.   rf\   r^    m    rri   vo    irt   ir\    ir^    *v^   ir»    /ts   ^"^    ^vi    /■~\   f^^    /^^   y\    ^-r    ir\  ^-^    i    tr\   r/^ 


o>Lrirn<rinrn  rooroco\Du-iinv£)MLncr<ocMor^o>0-3'io-vi-iHinoo 
lnoO'X)^-oooc^r-•oor^^^vo•<fu^u^o^r~o^roo10^r~o^r~•sl■cnoccvc 
o^-a•<■r-^-a■o^^AcN^nc;^~oo<roor^a^lnoo^r-c^lc^c^cN-^rvrooootn 
nc^«d■oc^^«d■'-^cooncrvo<^coroCQfor^l^lr^<^o<l•oc^Or^tH 
vocM<rvomc^i»iii— ^^■<r<^>3■l^<^s^«*lnc^^<t■<rr-^^o~a•<fifOl^u-l■<J•<■ 

^nli^u^u^L^lnl^,  lou^toi^'AiJ^inu^i/iinu",  ininiriinininioiT!  intnin 
I     I     I     I     I     t     I     I     I     I     I     I     I     I     I     I     I     t     I     I     I     I     I     I     I     I     I     I     I 

rHvoocrmcocN<rt^-v)-voincvic>-ioo\\ococ^jONLn(vioocNro<i-csiNvo 
o>omo-*cMvro>vfmcocrio-.  cNicv)i^ooovDoor-sr~c^iooto>3-a>tnvo 
tHrocr>vc>cviO'-H<ro>-<rr^mtnin<-HmocsioofOp^oooroi^rn>x)oc 
vOrH>-iinvDCOirir-ivocyi»d-o\rH<rr-»r<r*i-(tHr-iovoror-lconvooom 
vDcN  cn<J•vJ■«3■mvDvooc;o^c^ooor^cMLnmtnvo^.or~scoo^Or^^^c»^ 

C0C^C^l0^0^0^a^a^0^0^0^0^OOOOOOOOOOOOOIHl-^I-llH 

c^^cNCM(^Jcs^^Jc^^c^lcgc^lcNc^lMcn^^^n^n^o^Olrl^no^^n^n^^fn^^^^^o 
oocooocooooocooocooooooocooooooooooocooocooooooocooooooooo 


rH  ro  ro 
IT)  -^  ro 


CO  (N 


o  en  00 
r^  o  vo 


-*vooocr-~mo-icr><-occnvDLncnco 

C^Jr^COO^C<^Ot-^CMO^O>a•^-^t^vOOO 


r-.»3'cr>-<rrovocomsrcr>cMOfMcvJOcri<rr--<rK)-inovDr-(cncNooocM 

OOrHCr>OCNOrHi-lr-IOfnCNlCMrOCNO£NOi-HCNC^rHCJOOCOCr>rHfn 
^^  .-t.-Hi-lt-Ht-<f-(i-irHi-liHrHHrHi-liHi-HrH  i-Hi-1  i-HiH  iHiH 


(0 
0) 

O 


x> 


43 


CO         W   ti 


o 

S5 


c 

3 

to 


rHrocovorOrHo-sr-— icNcriu-icMvrc<nvocNjoo-3-iocNJcr>-vrcomocNionr~ 
vDCT\0'-<.-Ht-(cscMcnmu-)vDr~r^oooooc>)coco>crmr~coocnKr<i-vo 
iHiHcvicNicMcMcvjcsicNicNicvjcscsiesicNCMfnrofnforndmcn«3->a-sr'0-«3- 


69 


oooooooooooooooo 

o 

o 


^  vo  in  p~  "sT 
o  o  o  o  c 


vD  in  \0  -sf  \0  lO  vO 

o  o  o  o  o  o  o 


u 
o 
a. 

a 

(U 


<r  r~  CO  o 
o  m  -3-  cN 


evjin<)-ineNc>-ioocj>c>grHO»3-rHoeooi-Hr--mf^ooor~« 

I— IO^O^^DO^C^It'1«v^l^^Oln-a•C^^Of*^'^^OC^IfM^OCOO^P~• 


m  -a- 
rg  CM 


HooofOcoCT.cvirnaNfncom«3-v3-<)-<-<-o<r)rH<i-o«*-*<*ifnoooforo 


•o 


e^lr-^^oc^l^---cMlno^r-*^~-f>^a^^oa^l-^^^cNOf^oof^-<fOcotnoo^~c^lco 

r-l<rinOOCVJCNr-lrHi-I.Hi-ICN)cno4i-HiHCV|Vj-i-IOrHrHCSOrHOCMi-liH 

oooooooooooooooooooooooooooooo 
oooooooooooooooooooooooooooooo 


I  I  I  I  I  I  I  I 


I  I  I  I  I  I  r 


I 


I  I  I  I 


I  I  I  I 


CO 

o 

H 

o 

w 
o 


►J 
pa 

< 


(Li 

o 
!z: 

<J 

CO 

o 

H 
H 
w 
en 
o 

PLI 


ooOr^-*vOr^Olnfor~Ou^l-loom\DOOc^o^o^oo^o^oc^ll^Jc^Jc^^oo^ 

r^rOCOOOCOrHOOrHCv|iHOf^OfOCNaNOOCT>OOOCMCNiHONOOOrH 

ooOr-4t^Or^rHIHl^l-^r^r^t^t-^r^.HOr^r-^Or^OT^r^r^oOr^r^ 
oooooooooooooooooooooooooooooo 

I*   l'  I*  l'   I*   I*   I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I    ' 


rHOincoinr^-oorHr— ivocTii^o>c-jvoiofoooroc^iin«3'-<rc^cNjcg>3'cnr-i 
ONrOrHvoooiocsmcTi(Nir)mfoovooocOr-iONrHr~ooovDr~ooiHrHr^o 
o\vDr-or-~rH>*-<J-mocninHi-io>ootooovor-.cnt-Hoovoocvj-;i-r-~r~vo 
cv)f^eo<TirovOrHc^JvDor>-mr^oNrocnr~oooovo<NrHrHCN«3'0<TivDOoo 
ir>i-4vDr*mforofNja^r-~cvjvocvJOr-li-)i-irHf^)iHcocvic>'iroo>OrH<Tvi-iO 
^p^^cM>3•^^lln^a■<l•o•«d•■4■■*•<3•<^•a•-3■•3•«a•»3■^»a•«3•<rnfn^o-*^o 

ininioininininioiriu^ini/iinininininininioioirtiniriininininintn 
i.    I     I     i     I     I     I     I     I     I.    I     I     I     I     I     I     I     I     I     I     I     I     I     I     I     I     I     <     I     I 

v0O0^00O^CMv0Ou-|00r»lnc^^O^^0cn<N^OsI•00vO^~lHt^C^^~C^10^r^c0 

o\fnoo-*iAr~cM>x)r-iinrHrHcNi>a-coocviooom(NOi-Ho>ooo>mOLn 
Qr.»^^r)0^*1vD^^r^<l•^-ll/^o^»3■c3^-^•I-l^ocovDvoOr-lcoo^o^csvo^^r^ 
vOvOCMvOCr>OfOvOvOCTNOOOI~~r~OOOiHrH.HCv|rO-C-inininvOr^O>OiH 

coroioiniovor«.r-r~f-~ooooooo>a>ooooooooooooOt-ii-i 

^,Hr^IHI^■Hr^r^^^r^^^r^r^IHr-ICNC^I{^^C^)CVlCMC^JC^^t^JCMe^JC^^C^^C^JC^^ 

cocnrncnc'^MrncnfOfocncncocnf'ifOcocnrocnc^focncncnf'ifocoforn 
oooooooooooooooooooooooooooooooooooooooooooocooooooooooococo 


invocMr^^jeMvocMrnoocNJOvocM 
rH«d-<vjcvjcv)i-4aivoi-HCM-srunroo 


o> 

<■  O  C  rH 

o 

CM  00  -3-  fH  o 

cs 

00  CM  <r  o 

o 

m  o  o  «a-  c 

o 


rH    1-1    iH    CM   CM 


CM 


\r\  en  o  «* 

sraiinocoooco>3--*vosrr^r»«rHcoofor~-<rocMtHoOfOCMcvjr~mi-io 

oooi-Hi-ti^rorococMCMcMcMcocMmrHvocMoocnmocnocnvOCTicMco 


0) 
O 


CO 


.J 


in 

CO 

> 


p^   Ci   CO 
iJ  kJ  (J  hJ 


fd  <; 

m 

o 

pj  M  a 

H  cn 

— *   i-H 

rH 

»-H 

S  <!'-" 

a  s 

CD. CD 

CD 

CD 

o 

o 

z 

c 
a 


O 

vo  in  o 

invocMr-^oooiAvj-r^fHiOrHcorOpH^a-u-icMoocMvo 

vOvOO>CJNiHOCMCMrOCOcn-<r     I      1    OOOOOOOOOiCTi     I     ONOOOiHrHOCMCM 

•*<--;r«ji-LnmmminiriiriiriMwininu-immiocnmvo\o*i>vDvovDvov£> 


vovfor^OcvJcMvoo 
vovoo>CJNiHocMcMrococn-<r 


70 


oooooooooooooo 

o 

o 


ooooooooooooooo 


43 
O 

o 

w 

c 
cd 


00  iH  Oi  ro 

to  o  -d-  a^ 


rHinvoo\i— lo^e^^^-cnoool^lvDcno^^olA^-lsr^-^vDcoo 

000^-iriO<t'U^vDCT<-*<T>vDOrHr-ICTNOOOOOinf^vr-cr 


iH  00   Oi 

in  CO  CTi 


<rcn<rcMroc>inocno-o-rvJCT>cjNO>csc>iromcMOcncoocNCMOOrHCTi 

CTi- 


3- 


ON 


ooooMvo<rmocvicoo<Tivorooc>ioo«*c>jo><nrHoovof^'<tcvjoooo 

I^O•<^OrHOl^OOOOr^OfO^-lr^Ot-^l-^rHOCS^r^OI-^^^r-^CM!^<N 

OOOOOOOOOOOOOOOOOrHO  ooooooooooo 

oooooooooooooooooooooooooooooo 


I 


I    I 


I 


I    I    I    I 


I 


I 


I     I    I    r 


I    t     I    I     I     I 


in 


m 

z: 
o 

H 
H 
O 

w 

o 

w  pi 

tJ   Pm 

<  o 

H  23 

<; 

CO 

o 

H 
H 
H 

to 

O 
fl4 


M 
O 


ooooo^vo<^^^oop»I^l^o^occc^la^c^Jmmo^^^co<rcno•d■mo\Dvo•vt 

C;r^I-(^ncOr^OOC^OC^c^0^l-^>a■«3•rOC^lOr^O^OO^c^JOOO^CT^<fIH 

^^^^^r^r^H^^r^I^r^fHO.^.Hr^r^cMr^^^r-^^^o^-l.H.HOOr^H 
oooooooooooooooooooooooooooooo 


I    r    I     I     I     I    I    I     I     I     I    I    I     I     I    I    I    I     I'    I     I    < 


I     I    I    I    I    I     I 


<ror~o^Dr^iHcv|vo 

-a-r-.cMLnc>rHrofn<- 

iniHominv£)^-o\f*icr>r~» 

oo^rcNinomfOinoo 

or~»or-»oor^cyiOOcsro<Nr~ 

^fOfon(n'<tmrH>*csj<rfO 


00vOfM<■OO^O^^~v0m^0f1O00>^O^00C1I^O^C^^ 

roiH<rc^icNioo<>of-~for~c^ip«.ot^c^ioor-»"*'«o'n 

<ru-|^OOI^m<t^-^^^^^O^C^I«d■vOCnvOO^O^ 
r»-3-iHinmvfrHO<00'XlrHOCT>00C>>lCn00O 

>a-»*r-io>r~-coCTiin^i-imr~r~criiooooom 
<•lncN^r^^-^c^l-J•<r«a•c^J<^J-;)■LOcM•<t•lnr^m 


CO 

<r  oi  cs 


minmmmmmmminmininminiriininiotninLninmminminiom 
I     I     I     I     I     I     I     I     I     I     I     I     I     I     I     I     I    I     I     I     I     I     I     i     I     >     I     I     I     ' 

vor^ooin«3-ooin<r-*vocgmoofOvD-J-oor^cNooo>tHcocN)\or^-v»-oocoin 


vOm-J)-rHCOCTiCriCMrOOOOOrH<rO 

oo<Tim<y\vfvoinoorHcv)c^a>cr>cn 

I-^<MfO>3•a^O^I-lc^^rOfOvOr^OL^ 
i-li-<r-Hi-IOv|CVJCv|CNJCSrv4(nm 


vOvDOOCMr-^vO«3•r-~OOO^CT^^OCN^^^'^m 

iH<rvo-a-o>fnmcN)i^ovi-<-io<ricoCTi 
OrHcnst«d-miri>^>o-sfOmvovor~o 


(NC^«^lCNI<^JC^^^^JC^JC^^C^Je^l^^^C^lC^lCv^l^^C^^C^^C^lCVIC^^C^)(^JC^JC^IC^^CN<^JC^^C^^ 

cncocnfncocnrnro(^rocncoPicocnpicncoc*^fn(^rncnMrnoofncncnrn 
ocooooooooooooooooeooooooooooooooooooooooooooocooooooooooooo 


ooooo>^c4  0r^r--ooo^mlnmmlnlrl<rv£)OOlno^oo<r^o<^ln«3•^o 

O-d-  r~vOO  r~tHOOOOO<TiOCTi<T>OOtniHmvOCOOCMO<-<MCN 


<N 


P>1 


<0 


o  in  (^  00  tH  H 

oor~m<fincnr^oomCTNOocTsocN)CNjoot~~>«DP~fniHcn'nococN<roco 

(Nog(ncNC»JCO<yiCr>CNOCMCNlOCMfnrOCM.HiHCv4\OrOt-IC^CN)vOOOOfOH 


CO 
0) 

o 


> 


> 


00 

o- 

<  in 

^ 

fe^  pq 

CM      CO 

<CM 

CD    > 

m 

CD 

J3  VO 
Z  CO 


> 

2 


o 

55 

C 

U 

PQ 


CvJ 


r^oOrHcororHinfOoovocri>.ocviorovooofOcoin«a-oOrH-;f<j-c^or^f^co 
cNcM    I   co-3-inininininvot-^ooCTicTi<T>criOrHiHrHOr-(oo>v)-sr>>Dvoinr^ 


f-*i 


71 


minvominirimvDvoiOf~*minovD 

OOOOOOOOOOOOOrHO 

o 
o 


i-iOOOOOOOiH 


o 
o 

w 

c: 


vOCOmrO<x>e-l«*<-OfnvOaiCMU-|CSmiHf-^iHOCv|vDOOtH 

r^^oooI^o^r^ooo^vosrc^^c^^LnmI-lu^■<J•ooonooo^(^lHr^ 
oocomcMcNogcMfMocMcocMocroccnooo  o\oo  a\  oo 


p. 


vocMcoOr-<<tcMcMOcMrHmvou-ivoooooo\tNfnr~r^rxr^ 

00<-iHr-.u-|vOrHOiHOcO«NtHinvOiHO<l-r~.rHvDO\-d- 
OOOOOmoOOOOOOOt-Hr-IOOOOOOOOO 

oooooooooooooooooooooooo 


rill 


I    I    I    I    I    I    I    I    I    I    I    I 


o 

H 
H 

o 

OS 
»n  CM 

o 

M  Pi 
l-I    Ph 

« 

o 

o 


o 

p-< 


i-IO>0>OOCTiOO\00<-^CMCMrOCO>3-r-IOvOOOrHr-.ONCM 

c>i(NCNiiHvDoo<i-ooincri>a-coocvivoo<n<NmtHi-im-;fi-- 

Hi-lrHiHiriOlOrHiHOrHOt-ICMOOOOOCMiHOCMO 
OOOOtHOOOOOOOOOOOOOOOOOOO 


I       I      I       I 


I       I       I 


I 


I 


I       I      I       I       I 


I 


OiricNr-(fOcMroiri-3-cv)omo\u-icrimcri«3-cr>cncooOLnu-i 
vomp^<riOOvOt-iinoovDcri<yioovovoro-*cN<ro>f^in^.^ 

OOOmrHOmcv|VOOr-l<X)rHO\i-lsOeNr«»OOCMi-IOOCnr>. 
l^vOtHiHOCOoOcnvor^vDCTiu-ir^rHOOOOirivOrHrHncNvO 
mfn00Ov0vOv0i-IOOvDu-iO0>CTir~>r0cSOvDvO<J-iHO 

rHcvj^u-)cof>.jinfncvjinrHvo.<rvoovo-sj-«*mFHm.Hincvj 

I     I     I     I     I     I     I     I     I     I     I     I     I     I     I     I     I     I    I     I     I    (     I     I 

a>r-iH<yiCTiOincNrHr^or~t-icv)CT>cr>ooornr-ivD»nLno 

f^OOvOa3CfncN|CNJCMinrOrHCOCOP^rHrHair~-;)-CNCMvOCvJ 

r^iri«;fv£><)-o»3-oovrocvjocv)cMinii-ioN<roo«a-rvivrvi30 
voc»i»*oro<rrHcscovoo\cMcriroa^a\ooCTi<rtncT\irit-»cM 
i^ooOiHrHMeMr^co<)-vd-vDr--oocN)<i-cr>iH.3-^<t-\ooor^ 

CO0O0OO00OCOCX30O0O0O0O000O0OOO0OCO0OOO0OOO0OO00O 


cNai>£5roocvi<-t>«ooo^rHiHfncNvoror^<}-oooocN|[no 
OOoo-Nrr^i^CTiiTio>r~^iHOOrH-d-r-<o<j-vooMcM«d-r-r^ 


00 

ca 


I    I 


I    I 


HO  or^a>o-3-  s;t 

MavncnrH^ooinrHooinrHoocMmrHcNOiHcMo«^o 

o^o^c^I^r^r^c^)c^4c^loc^J^^a^^^a^voo^r-^oo^Or-^oH 


(0 
O 


<!   nj   n)  to 


o 
en 

CO 


•a- 
en  CO  m 

ro  fn  CO 
rt  ro  1-1  CO  ^ 

(0 


00  r^  rt  (0 

CO  CO  CO  <r 

CO  CO  CO  CO 

iH  iH  cy\  a> 


o 

« 

a 

9 
M 
pq 


<r  00  vo  ON  lo  CM  r^ 

vOvOOiHCOCMOOO«*<J-rs.l^OOOOOOCMOi-<Ht-ICMCOin 

f^oocMcocgco<t-i^vx>ooooocMcsmoocj>ooooooo 
r^i^oooocoooooooooooooa»<r>cy.cr.cr>aiiHrH>-ii-ii-ir-i,-4 


72 


n 

u 

•r4 

rt    1    to 

M 

0) 

ti 

0)   -H 

CO 

•3 

O 

>>    V4     CO 

VI 

tI 

•H 

^   o   rt 

X   rt 

CO 

«J 

M-l  43 

o  -u 

(3 

Id 

t3 

Pn     to 

O 

«i 

0)     Q)    01 

O 

o 

W   H  J= 

•3    3 

a 

rt  ^  +J 
C   rt 

3    <U 

rt  M-i 

a 

60^    (3 
•H    O    O 

3    rt 

3. 

•0 

(0    »4 

o 

O 

d 

<U    04  U 

.H    U 

Rt 

•TJ        rt 

r-i     O 

• 

73 

u 

rt  4J 

rt  4-1 

o 

3 

4J 

u        to 

b 

• 

rt 

CO 

<U    (0 

#» 

-* 

4J   <u  -d 

>»   M 

1^ 

• 

/-N 

c 

4J    4J    c 

J3     O 

CTi 

X 

to 

t( 

rt   o  3 

rH     C    O 

>^^ 

rH 

3. 

u 

tU    M 

n  r«» 

J3 

u 

to 

0) 

a)   no    00 

4J  m 

O 

00 

,a 

J3           (U 

Q)  0> 

O 

o 

»« 

a 

w   rt   M 

e  -H 

P. 

jj 

4J 

a 

N^                     O 

o  ^ 

(U 

o 

X 

c 

IM 

u 

J3 

01 

M      • 

O    3 

ft 

&. 

u 

m 

<U  /-s    P^ 

ji  o 

o 

^   "V   r-i 

(X  (0 

m 

UH 

3 

01 

B     O  XI 

3 

en 

0 

•H 

J3 

3   +J   rt 

0)  j: 

rH 

4J 

C  -u  Ja 

*J    o 

to 

-3 

•H    O 

o  *-} 

3 

<U 

0) 

A 

a  6  u 

C 

o 

•H 

3 

M 

pa   o  Ou 

(U    B 

•H 

M 

•rl 

0) 

T>    O 

4J 

<U 

rt 

•O 

M  =     rt 

M 

rt 

(0 

H 

s 

o    • 

to  iw 

4J 

CU 

9 

0     to 

J2 

3 

4J 

X 

a 

«m  (U 

0)      1      4-1 

4J       - 

3   m 

(U 
•H 

3 
(U 

<u 

a 

B  =     o 

<U  ,3 

M 

O 

to 

9 

rt        G 

• 

u  u 

o 

0) 

rt 

M 

a  to  Hi 

/-N 

-3 

M 

« 

C  -3  CT> 

3    (U 

u 

3. 

(U   -H 

\o 

•H     »4 

o 

0) 

d 

fH   'O  J3 

a^ 

-3 

UH 

JS 

(U 

^    (U 

rH 

CO    3 

+J 

,3 

<M 

rt   u 

(U    3 

3 

4J 

o 

•H    (I)      • 

M 

•O  ^ 

o 

«+H 

(-1    M    3. 

(U 

3 

•H 

o 

MH 

lU 

rt    Cu 

A! 

4J    3 

4-> 

O 

o 

>        iw 

H 

•H   -H 

rt 

J3 

a 

0)    O 

rt 

3 

3 

o 

u 

01 

"rC 

S 

oo  CO 

•H 

o 

o 

<a 

M    -U    CO 

^-x 

(0  tu 

H 

&. 

u 

x> 

0)          -H 

e  -3 

O 

(U 

u 

CO 

• 

4J  ^    to 

U 

3 

(U 

<u 

5 

.U   4J    rt 

o 

4J 

•3 

3 

(U 

W 

0)   'iH   XI  rH 

••   -H 

rt 

3 

J= 

•H     S 

o 

M    3 

• 

-3 

(U 

a 

4J 

QJ 

u 

O    60 

/-s 

3 

B 

(U 

>> 

j«!  M  x: 

rH    rt 

o> 

rt 

B 

Ci 

.D 

(U    0)    *J  "3 

o   S 

vO 

(U 

H 

0)   ^ 

3 

a 

o> 

3 

x: 

H 

•o 

M    6    3 

rt 

.» 

H 

0 

u 

rt 

01 

O    3    O 

X)  /-^ 

V— ' 

•H 

B 

•  • 

•o 

3 

(U 

C  -a 

0) 

M 

M 

v 

<U          H  T) 

a    (U 

M 

3 

•  • 

O 

(U 

o 

xs  -u  rt 

3 

j: 

(U 

<U 

x: 

MH 

^ 

(U 

H   -H    *J 

4J 

<U    CO 

^ 

a 

u 

0 

M 

bO  CO 

•H 

T3   -H 

H 

tn 

o 

0) 

3 

D. 

•H 

3 

3    rH 

rt 

rt 

(X 

x: 

Z 

••  "3  "3 

60 

+J    X> 

S 

(U 

w 

a 

M           C 

rt 

•H    3 

u 

• 

C 

o 

tU    h    3 

0 

3  a  e 

X 

3 

rH 

3 

> 

4J    3    O 

bO  3 

o 

60 

rt 

•  •        rt 

U 

•H 

O     O    H 

i|^ 

rt    3 

u 

•H 

0) 

PL       3 

PQ 

00 

Z   iw    60 

o 

S  •-' 

UH 

OS 

S 

O         a" 

.a 

u 

•3 

(U 

UH               OJ 

APPENDIX: 
DETAILED  SUMMARY  OF  PLATE  DATA 


73 


74 


TABLE  A-1 

LIST  OF  PLATES  USED  IN  THIS  INVESTIGATION 

(Numbers  are  in  the  system  of  the  respective  observatories) 


It 


Allegheny 


15  572  (1919)  103  394  (1969) 

15  573     ••  103  512  (1970) 

15  684     "  103  522 

15  685     "  104  157 

28  668  (1922)  104  295  (1971) 

62  816 


McCormick 


17  148  (1924)                        77  407  (all  epochs  1968) 
196  "  ^08 

228  "  580 

2  29  "  5  81 

865  (1925)  624 

866  "  625 

966  "  650 

967  "  651 

18  905  "  698 
906  "  699 
92  2  "  7  26 
923  "  727 

19  480  (1926)  744 
481  "  745 

564  "  812 

565  "  813 

20  573  "  814 
742  "  815 
74  3  "  844 
82  9  "  84  5 

21237"  866 

238  "  867 


Y  e  r  k  e  s 

0    4    (1905) 
0    7     (1905) 

012  (1905) 

013  (1905) 

0    18    (1906) 
0    30    (1907) 
0    52    (1908) 
0432    (1922) 

Z6-Se 

"Neb.     d'Orion' 
1    (1902.' 

iia 
913) 

4    (1914.120) 

-1^^ 
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TABLE    A-l--contlnued 


Zo-Se--continued 


2  (1910.110) 

3  (1911.056) 


5  (1916.069) 


San  Fernando 


1490  (1893) 
1496  (1893) 
3890  (1906) 


University  of  South  Florida 


13^  (all  epochs  1974) 

15 

16b 

22 

28^ 

29b 

44 


45 
55 

57b 

60 

63b 

65b 


Sproul 


T346b 


^The  identification  given  in  Chevalier  (1933)  .   The 
times  are  given  to  three  places  to  aid  in  identification. 

b. 


Treated  as  two  plates. 
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No, 


TABLE  A-2 
DATA  FOR  USF  OBSERVATORY  PLATES 


Serial  a 


Date         Grating    Exp.    Center       HA 


13^ 

9 

Mar 

1974.189 

- 

7 
2 

min  . 

1 

+ 

3  min. 
5 

15 

9 

Mar 

1974.189 

- 

2 

1 

+  156 

16^ 

10 

Mar 

1974.192 

— 

12 
2 

1 

+  '  ' 

+ 

6b 
67 

22 

U 

Mar 

1974.203 

- 

9 

1 

+  103 

283 

24 

Mar 

1974.230 

2 

9 
1. 

5 

1 

+  120 

+  125 

29^ 

24 

Mar 

1974.230 

- 

9 

1. 

7 

1 

+  140 
+  146 

44 

17 

Sep 

1974.710 

- 

9 

1 

- 

93 

45 

29 

Sep 

1974.743 

- 

8 

1 

— 

35 

55 

17 

Nov 

1974.876 

"■ 

1. 
4. 

5 
5 

2 

— 

6 
3 

57^ 

17 

Nov 

1974.876 

— 

5 
1 

3 

+ 

+ 

60 
66 

60 

18 

Nov 

1974.879 

— 

1 
6 

2 

: 

12 

7 

63^ 
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Opiate  centers:   la   =  5  h.  48  min.,  6^  =  -5°,  24' 

(Trapazium) 

2  30'  N,  1  min.  W  of  Trapezium 

3  30'  S,  1  min.  E  of  Trapezium 

''separate  exposures  reduced  as  separate  plates. 
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